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White  noise  spectra  of  diodes  1N746  to  1N752  have 

been  used  to  investigate  the  reverse  conduction  mechanisms 

in  silicon  junctions  breaking  down  by  a  combination  of 

tunneling  and  impact  ionization  multiplication.   The  reverse 

voltages  at  the  onset  of  multiplication  noise  indicate 

that,  regardless  of  junction  field  strengths,  carrier 

multiplication  sets  in  only  when  the  carrier  energies 

exceed  the  threshold  energies  for  ionization.   Further 

examination  of  noise  data  yields  their  values,  indicating 

that  the  multiplication  process  is  dominated  by  the 

effects  of  the  threshold  energies  for  ionization.   When 

they  are  taken  into  account,  it  becomes  clear  why  Tager's 


and  Mclntyre's  theories  of  noise  in  avalanche  diodes, 
and  the  conventional  DC  multiplication  calculations,  are 
not  applicable  to  low  breakdown  voltage  diodes.   Con- 
sequently, a  new  method  of  performing  DC  multiplication 
calculations  as  well  as  a  new  theory  of  noise  in  low  break- 
down voltage  diodes  is  developed. 


CHAPTER  I 
INTRODUCTION 

As  the  reverse  bias  across  a  p-n  junction  is  increased, 
the  junction  eventually  begins  to  pass  appreciable  current. 
Typical  voltage-current  characteristics  for  commercially 
available  diodes  are  shown  in  Figure  1.   It  is  evident  that 
the  reverse  current  characteristics  are  considerably  dif- 
ferent for  p-n  junctions  which  "break  down"  at  different 
voltages. 

The  type  of  breakdown  typified  by  Figure  la  is 
characterized  by  (i)  temperature  sensitive  pre-breakdown 
(saturation)  current  which  approximately  doubles  for  each 
10°C  rise  in  temperature,  (ii)  a  well-defined  breakdown 
voltage,  V  ,  and  (iii)  a  unique,  positive,  breakdown 
voltage  temperature  coefficient,  B,  defined  by 

VB  =  VBOf1+B(T-T0)] 
where  VBQ  is  the  breakdown  voltage  at  room  temperature,  T  . 
This  type  of  breakdown,  characteristic  of  relatively  wide 
p-n  junctions  with  V  >14  volts,  is  commonly  known  as 
"avalanche"  breakdown  since  it  can  be  understood  in  terms 
of  an  impact  ionization  mechanism  resembling  an  avalanche 
(Tyagi,  1968)  . 


Reverse  conduction  of  the  type  shown  in  Figure  lc 
is  characterized  by  (i)  temperature  insensitive  reverse 
current  —  temperature  change  of  200°C  changes  the  current 
by  a  factor  of  two  or  three,  (ii)  lack  of  well-defined 
breakdown  voltage,  and  (iii)  a  negative  temperature 
coefficient  whose  magnitude  depends  on  temperature  and 
reverse  voltage.   This  type  of  breakdown,  characteristic 
of  narrow  p-n  junctions  capable  of  passing  considerable 
current  densities  at  reverse  voltages  smaller  than  about 
2.5  volts,  can  be  understood  in  terms  of  internal  field 
emission  of  valence  band  electrons,  also  known  as  Zener 
tunneling  (Chynoweth   et  al. ,  1960;  Tyagi,  1968). 

P-n  junctions  Wiiich  break  down  between  about  2.5  to 
14  volts  do  so  as  a  result  of  the  combined  effect  of  internal 
field  emission  and  impact  ionization  multiplication  (Tyagi, 
1968) .   Since  most  of  the  measurements  which  formed  the 
basis  for  the  understanding  of  p-n  junction  breakdown  were 
DC  measurements,  this  range  of  breakdown  voltages  has 
never  been  carefully  investigated.   The  reason  for  it 
was  the  severe  difficulty  encountered  in  distinguishing 
the  contributions  of  each  of  the  two  mechanisms. 

The  aim  of  this  study  was  to  investigate  the  reverse 
conduction  mechanisms  in  diodes  breaking  down  in  this 
range  of  voltages.   The  difficulty  which  previous  workers 
encountered  in  distinguishing  the  tunneling  current  from 
the  impact  ionization  multiplied  current  was  overcome 


by  using  noise  as  a  diagnostic  tool,  the  idea  being  that 
noise  would  more  closely  reflect  the  basic  microscopic 
conduction  mechanisms.   This  technique  was  first  employed 
by  Lauritzen  (1966)  in  determining  the  onset  of  carrier 
multiplication  in  low  breakdown  voltage  silicon  diodes. 
Interest  in  avalanche  diodes  as  photodetectors ,  micro- 
wave signal  sources,  and  amplifiers  has  already  led  to 
extensive  experimental  (Haitz,  1965,  1966;  Haitz  and 
Voltmer,  1968;  Baertsch,  1966,  1967;  Kumo   et  al.,  1968) 
and  theoretical  (Tager,  1965;  Mclntyre,  1966;  Hines,  1966; 
Gummel  and  Blue,  1967)  studies  of  noise  in  these  diodes 
at  low  and  high  frequencies. 

Because  noise  measurements  provide  the  experimental 
foundation  of  this  study,  Chapter  II  is  primarily  devoted 
to  the  experimental  considerations  and  techniques  involved 
in  collecting  the  noise  data.   DC  and  capacitance  measure- 
ments are  discussed  there  also. 

A  general  discussion  of  the  data  is  undertaken  in 
Chapter  III.   It  is  found  that  Tager' s  (1965)  and  Mclntyre' s 
(1966)  theories  of  noise  in  avalanche  diodes  are  incapable 
of  explaining  the  data.   The  reason  for  this  becomes  clear 
when  the  mechanism  of  the  onset  of  carrier  multiplication 
is  inferred  from  the  noise  data.   It  is  then  also  apparent 
that  the  conventional  method  of  calculating  DC  multipli- 
cation is  not  applicable  to  low  breakdown  voltage  diodes. 


When  a  new  definition  of  the  DC  multiplication 
factor,  and  a  new  functional  form  of  the  expression  for 
noise  are  introduced  in  Chapter  IV,  the  noise  data 
reveal  in  detail  the  influence  of  the  threshold  energies 
for  ionization  on  the  multiplication  process.   On  basis 
of  the  inferences  made  in  this  and  the  preceding  chapter, 
new  theories  of  noise  in  low  breakdown  voltage  diodes 
are  developed  and  compared  with  experimental  results. 

The  details  of  the  multiplication  process  inferred 
in  Chapters  III  and  IV  are  incorporated  in  Chapter  V, 
which  discusses  a  simulation  program  implementing  a  new 
method  of  performing  DC  multiplication  calculations  in 
low  breakdown  voltage  junctions.   The  results  of  the  con- 
ventional and  the  new  calculations  are  then  compared  with 
data. 

Chapter  VI  summarizes  the  highlights. 

Chapter  VII  suggests  topics  for  further  study. 


CHAPTER  II 
EXPERIMENTAL  METHODS 

The  p-n  junctions  investigated  in  this  study  are  com- 
mercially available  diodes,  types  1N746  to  1N752,  manufactured 
by  alloying  small  aluminum  buttons  onto  n-type  silicon.   The 
raised  appearance  of  the  aluminum  buttons,  shown  in  a  vertical 
cross-section  of  the  diode  in  Figure  2,  is  the  result  of  acid 
etch  used  to  delineate  the  junction  geometry  and  eliminate 
surface  effects  around  the  button  periphery. 

The  samples  chosen  for  study  were  those  which  exhibited 
very  little  flicker  noise.   In  final  samples,  flicker  noise 
was  not  noticeable  at  frequencies  greater  than  300  Hz.   This 
selection  criterion  presumably  eliminated  from  consideration 
those  diodes  in  which  surface  or  edge  effects  might  have  had 
significant  influence  on  the  reverse  conduction  mechanisms 
of  the  junctions.   All  final  samples,  whose  DC  V-I  charac- 
teristics are  shown  in  Figure  3,  came  from  Continental  Device 
Corporation,  now  Teledyne  Semiconductor. 

In  this  study  the  reverse  conduction  mechanisms  were 
inferred  from  two  sets  of  diode  terminal  measurements:   the 
DC  V-I  characteristics  and  noise.   The  DC  characteristics 
provided  information  on  the  average  behavior  of  carriers  in 
the  junctions,  whereas  the  noise  measurements  provided 


statistical  information  on  the  transient,  or  microscopic, 
details  of  the  conduction  process.   The  noise  measurements 
thus  provided  the  groundwork  details  for  the  conduction 
model,  which  will  be  discussed  in  the  following  chapter. 

The  noise  measurements,  which  form  the  experimental 
foundation  of  this  study,  were  performed  with  the  aid  of  the 
noise  measurement  system  shown  in  Figure  4.   The  GR  1381 
Random  Noise  Generator  coupled  to  the  HP  350-D  600  Q   attenuator 
supplied  a  2  Hz  to  50  KHz  white  noise  calibration  signal. 
The  600  n  resistor  following  the  attenuator  provided  the 
proper  impedance  match  for  the  attenuator,  while  the  1  pf 
capacitor  in  series  with  the  200  Kfl  resistor  provided  DC 
and  impedance  isolation  from  the  rest  of  the  circuit.   The 
impedance  isolation  was  necessary  to  maintain  a  constant 
impedance  level  at  the  attenuator  terminals  regardless  of 
the  impedance  levels  existing  in  the  diode  bias  network, 
and  to  convert  the  noise  calibration  network  into  a  high 
impedance,  current-like   source  to  prevent  excessive  loading 
of  the  diode.   The  diode  bias  network  consisted  of  a  large 
number  of  wire  wound  resistors,  R,  .   When  used  to  adjust  the 
diode  reverse  current,  they  were  chosen  so  that  R  was  more 
than  ten  times  greater  than  the  diode  dynamic  resistance,  R , . 

The  low  noise  preamplifier  used  for  its  first  stage  a 
preselected  sample  of  a  Texas  Instruments  low  noise  JFET ,  the 
SFB  8558.   The  schematic  diagram  of  the  preamplifier  is  given 
in  Figure  5,  and  the  preamplif ier ' s  equivalent  noise  resistance 


referred  to  the  input,  R  ,  as  a  function  of  frequency,  is 
given  in  Figure  6.   The  preamplifier  was  followed  by  a 
selectable  gain  amplifier  which  amplified  the  noise  signal 
to  a  level  suitable  for  processing  by  the  GR  1925-1926 
Real  Time  Spectrum  Analyzer.   The  oscilloscope  was  used  for 
visual  observation  of  the  noise  waveforms. 

The  GR  1925-1926  Real  Time  Spectrum  Analyzer  contains 
a  bank  of  45  third-octave  filters,  ranging  in  center  frequencies 
from  3.15  Hz  to  80  KHz.   The  output  of  each  filter  is  sampled 
for  up  to  32  seconds  and  the  dB  of  the  RMS  voltage  of  each 
filter  is  computed  and  displayed  on  the  GR  1926,  or  is  printed 
out  by  the  MDS  800  tape  printer. 

The  measurements  were  performed  as  follows.   First,  the 
noise  calibration  signal,  provided  by  the  GR  1381  and  the 
600  ft  attenuator,  was  removed  by  disconnecting  the  attenuator 
from  the  circuit  and  replacing  it  with  a  600  ft  resistor 
placed  in  parallel  with  the  600  ft  attenuator  termination 
resistor.   The  diode  bias  resistor,  R,,  was  then  chosen  to 
obtain  the  desired  value  of  diode  reverse  current,  and  five, 
32  second,  diode  noise  data  sets  were  recorded.   Next,  the 
600  ft  resistor  in  parallel  with  the  attenuator  termination 
resistor  was  removed,  and  the  attenuator  returned  to  the 
circuit.   The  attenuation  level  was  then  adjusted  to  obtain 
a  calibration  plus  diode  noise  output  about  20  dB  higher 
than  the  diode  noise  output  alone,  and  a  second  set  of  five, 
32  second,  readings  was  recorded.   Given  the  system  noise 
equivalent  input  circuit  of  Figure  7,  it  i s  a  trivial  matter 


to  determine  that  the  diode  noise  current  spectral  density, 
S.,,  is  then  given  by 
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where  S    is  the  spectral  density  of  the  noise  generator  out- 
put voltage,  R   is  the  equivalent  noise  resistance  of  the  low 
noise  preamplifier,  R   is  the  diode  dynamic  resistance,  R, 
is  the  bias  resistance,  C   is  the  total  input  capacitance, 
dB     is  the  attenuator  setting,  and  AdB  is  the  difference 
in  dB  between  the  calibration  plus  diode  noise  level,  and 
the  diode  noise  level  alone.   This  calculation  was  carried 
out  for  each  filter,  and  the  final  results  were  averaged  over 
all  filters  processing  a  white  noise  spectrum.   Typically, 


this  resulted  in  averaging  over  the  25  filters  of  center 
frequencies  greater  than  250  Hz. 

The  system  accuracy  was  verified  by  measuring  the  thermal 
noise  of  resistors.   When  resistors  of  200  fi  to  2  Ml]  were 
inserted  in  place  of  a  diode,  and  the  above  measurement  pro- 
cedure was  repeated,  the  resistance  values  predicted  from 
noise  measurements  agreed  to  better  than  2%  with  values 
obtained  from  precision  bridge  measurements. 

Since  the  low  noise  preamplifier  employed  a  JFET  for 
the  first  stage,  it  sensed  the  voltage,  rather  than  the  current, 
at  its  input  terminals.   Consequently,  the  diode  noise  current 
was  converted  to  a  noise  voltage  by  the  diode  dynamic  impedance 
before  the  noise  signal  was  amplified.   The  diode  dynamic 
impedance,  therefore,  had  to  be  determined,  and  was  measured 
using  the  diode  bias  network,  capacitively  coupled  to  a  Wayne 
Kerr  B601  Radio  Frequency  Bridge  utilizing  a  Wayne  Kerr 
SR  268  Source  &  Detector,  as  shown  in  Figure  8.   The  100  KHz 
signal  applied  by  the  bridge  to  the  diode  was  adjusted  until 
further  signal  level  reduction  produced  no  difference  in  the 
measured  RC  values.   This  typically  occurred  for  applied 
signal  levels  of  less  than  4  0  mV  RMS.   The  parallel  equivalent 
RC  values  of  the  diode  and  its  bias  network  were  thus  deter- 
mined to  an  accuracy  of  ±1%. 

The  capacitance  values  determined  above  were  also 
employed  in  the  C-V  plots  used  in  estimating  the  junction 
doping  profiles  and  doping  densities.   The  diode  areas  wore 
obtained  from  junction  photographs  taken  after  the  aluminum 


button  was  etched  away  in  phosphoric  acid.   The  error  in 
determining  the  junction  area  in  this  manner  was  estimated 
to  be  about  15%. 

A  typical  junction  surface  is  shown  in  Figure  9.   Although 
the  surface  is  full  of  ridges,  it  is  assumed  that  the  junc- 
tions are  very  nearly  planar  because  they  are  extremely  thin. 
On  the  scale  of  the  several  hundred  angstrom  thin  space 
charge  regions  (SCR) ,  the  ridge  contours  should  appear  as 
gentle  undulations  in  the  junction  topology. 

The  DC  V-I  measurements  were  made  using  the  test  setup 
shown  in  Figure  10.   It  consisted  of  the  diode  bias  network, 
a  Keithley  Model  615  Digital  Electrometer  operating  as  an 
ammeter,  and  a  Fluke  Model  801B  Differential  Voltmeter.   The 

accuracy  of  the  DC  measurements  was  ±0.5%  for  currents  greater 

-7  -7        -9 

than  10   amperes,  and  ±2%  for  currents  between  10   and  10 

amperes.   Voltage  measurements  were  accurate  to  +0.05%. 

All  measurements  were  taken  at  room  temperature,  22±2°C. 


CHAPTER  III 
DISCUSSION  OF  EXPERIMENTAL  RESULTS 

In  order  to  make  the  forthcoming  discussions  of  noise 
and  DC  V-I  measurements  meaningful,  it  is  necessary  to 
know  the  junction  properties  and  parameters.   Some  of  them 
can  be  inferred  from  the  -r  vs.  V   plots  for  the  four 
most  important  diodes,  the  1N749  to  1N752,  shown  in  Figures 
11  to  14. 

Since  the  — ~  vs.  V  plots  do  not  yield  straight  lines 
over  the  entire  range  of  reverse  voltages,  it  appears  that 
the  junctions  are  not  abrupt.   Unfortunately,  —~   vs.  V 

c      r 

plots  do  not  yield  straight  lines,  either.   Consequently, 

it  seems  that  the  junction  profiles  are  composite  of 

linearly  graded  and  abrupt.   If  it  is  assumed,  however, 

that   to  a  first  order  approximation   the  diodes  can  be 

treated  as  step  junctions  (provided  that  the  doping  densities 

are  evaluated  from  — ^  vs.  V   slopes  obtained  for  large 

CZ  r 

reverse  biases,  where  the  doping  densities  are  most  likely 

to  be  representative  of  the  bulk  doping  densities)  the 

results  of  Table  I  are  obtained. 

The  calculated  doping  densities  on  the  n-side  of 

the  junctions  are  seen  to  lie  within  the  range  of  doping 

densities  predicted  from  the  starting  resistivities  of 
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the  n-type  wafers.   This  seems  to  indicate  that  for  reverse 

voltages  greater  than  about  1.5  to  2  volts,  the  doping 

densities  at  the  junction  edges  attain  the  bulk  doping 

densities,  and  the  junctions  should  begin  to  assume  step 

junction  characteristics.   The  changes  in  the  — ^  vs.  V 

C 
slopes  for  smaller  reverse  biases  are  caused  by  the  extreme 

thinness  of  the  diode  SCR's:   as  the  reverse  biases 

decrease,  the  SCR's  narrow  and  enter  the  diode  metallurgical 

transition  regions  where  the  doping  densities  are  no  longer 

constant. 

However,  demonstrating  that  the  doping  densities  at 
junction  edges  reach  bulk  doping  densities  does  not  con- 
clusively prove  that  for  V  greater  than  1.5  to  2  volts 
the  junctions  may  be  approximated  by  step  junctions. 
Another  piece  of  evidence  which  strengthens  this  contention, 
though,  will  be  obtained  when  the  DC  V-I  characteristics 
are  examined.   This  discussion  will  be  most  meaningful 
after  the  noise  measurements  are  examined. 

The  typical  noise  data  obtained  for  diodes  1N746 
to  1N752  are  shown  in  Figures  15  to  21.   The  diode  noise 
current  spectral  density,  S. ,,  has  been  normalized  with 
respect  to  the  noise  current  spectral  density  of  a  saturated 
thermionic  diode  conducting  a  DC  current  equal  to  the  p-n 
junction  reverse  current.   The  noise  ratio,  MR,  is,  therefore, 
given  by 


NR  =  ^if-  (3.1) 

where  I   is  the  p-n  junction  reverse  current. 

Since  the  noise  current  of  a  saturated  thermionic 
diode  arises  from  the  reception  at  the  anode  of  all 
individually,  and  randomly  emitted,  thermally  excited 
cathode  electrons,  a  noise  ratio  of  unity  indicates  that 
the  p-n  diode  reverse  current  is  also  due  to  collection 
of  randomly  and  singly  generated  carriers,  or  carrier 
pairs.   This  is  precisely  what  happens  when  the  p-n 
junction  reverse  current  arises  from  internal  field 
emission  in  which  thermally  excited  valence  band  electrons 
tunnel  to  the  conduction  band.   The  noise  measurements 
thus  furnish  a  foolproof   method  of  determining  the  range 
of  reverse  voltages  for  which  the  diode  reverse  conduction 
arises  from  a  single  step  tunneling  process.    (A  few 
samples  exhibiting  multi-step  tunneling  were  also  observed. 
They  are  characterized  by  a  noise  ratio  less  than  unity 
(van  der  Ziel,  personal  communication),  as  shown  in  Figures 
21a  and  21b.   These  samples,  however,  exhibited  substantially 
greater  —  noise  than  the  remaining  samples  and,  therefore, 
were  excluded  from  further  study.) 

An  independent  method  of  verifying  that  over  the 
range  of  reverse  voltages  for  which  the  noise  ratio  is 
unity  the  reverse  current  arises  from  single  step  tunneling 
is  obtained  from  the  laws  of  internal  field  emission. 
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Chynoweth   et  al.  (1960)  and  Tyagi  (1968)  found  that  for 
a  phonon-assisted  tunneling  process  the  diode  reverse 
current  is  given  by 


I   =  AV"E  e  ;  Y  =  —  (3.2) 

r      r 

3qh 

where  I   is  the  diode  reverse  current,  V   is  the  diode 
r  r 

reverse  voltage,  A  is  a  constant  for  a  given  temperature, 

E  is  the  electric  field,  E   is  the  band  gap  energy,  m 

is  the  effective  mass  of  tunneling  electrons,  q  is  the 

electronic  charge,  n=-=— ,  where  h  is  the  Planck  constant, 

YE  3/2 
-(-I— ) 

and  n  and  n  are  constants.   The  factor  e  represents 

the  phonon-assisted  tunneling  probability,  E   accounts 
for  image  fields  or  Coulombic  forces  corrections  included 
in  the  theory,  and  V   takes  into  account,  in  an  empirical 
way,  the  effect  of  electron  tunneling  from  the  conduction 
band  back  to  the  valence  band  so  that  for  zero  reverse 
bias  the  net  current,  I  ,  is  also  zero.   For  abrupt  junction; 
Tyagi  claims  n-1,  n-1.5. 

YE  3/2 

(    g      ) 

Compared  to  e         ,  AV  E   is  a  slow  varying 
function  of  the  reverse  voltage  so  that  the  reverse  current 
may  be  written  as 


(3.3) 
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where  I   and  s   are  considered  constant  provided  V>0, 
o  r 

and  E  is  an  "effective"  electric  field  such  that 


s     w    s 
e_r"  =  /   e_FT3n"dx  (3.4) 

o 

where  the  integration  extends  over  the  entire  width  of 

the  diode  space  charge  region,  w.   For  an  abrupt  junction 

the  main  contribution  to  the  integral  comes  from  electric 

field  values  close  to  the  maximum  electric  field  value,  E^. 

Thus,  little  error  is  made  in  assuming  that  the  "effective" 

electric  field  is  given  by  E  .   Consequently,  for  an  abrupt 

junction 


/v* 


(3.5) 


where  s  is  a  constant  and  <J>  is  the  built-in  potential 
of  the  junction.   Thus,  if  a  junction  is   indeed   abrupt 
and  the  reverse  current  is  due  to  tunneling,  log  Ir  vs. 


should  give  a  straight  line. 

Inspection  of  Figures  22  to  28,  which  show  log  Ir 
vs.    1    for  diodes  1N746  to  1N752,  reveals  that  over 


r 


ft 

the  range  of  reverse  voltages  for  which  the  noise  ratio  is  unity 
log  I   vs.    1         is,  indeed,  a  straight  line.   The  DC 


V-I  data  and  the  noise  data  are,  therefore,  consistent, 
and  the  contention  that  the  junctions  may,  to  a  first  order 
approximation,  be  considered  abrupt  is  further  strengthened. 
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If  the  randomly  emitted  tunneling  electrons  undergo 
multiplication  due  to  impact  ionization,  the  randomness  of 
the  multiplication  process  generates  additional  noise, 
and  the  noise  ratio   becomes  greater  than  unity.   Comparison 
of  Figures  15  to  21  with  Figures  22  to  28  reveals  that 
over  the  range  of  voltages  for  which  the  noise  ratio 
increases  from  unity,  log  I   increasingly  departs  from 
the  straight  line  it  followed  at  lower  reverse  voltages. 
This  suggests  that  the  failure  of  the  reverse  current 
to  obey  the  tunneling  relation  is  due  to  impact  ionization 
multiplication  of  the  tunneling  current. 

If  this  contention  is   indeed   correct,  then  it 
should  be  possible  to  verify  it  quantitatively  with  the 
help  of  noise  theories  of  Tager  (1965)  or  Mclntyre  (1966) 
relating  the  spectral  density  of  noise  current  generated 
in  impact  ionization  multiplication  to  the  terminal  DC 
current  flowing  through  the  diode.   According  to  Tager, 
the  low  frequency  spectral  density  of  noise  current 
produced  in  an  impact  ionization  multiplication  process 
in  which  the  electrons  and  holes  have  the  same  ionization 
coefficients,  is  given  by 

Si  =  2qIt(M')3  (3.6) 

where  I  is  the  current  initiating  the  multiplication 
process,  which  in  this  case  is  the  tunneling  current, 
and  M   is  the  DC  current  multiplication  factor,  defined  by 


(3.7) 


where  I   is  the  terminal  reverse  current  flowing  through 

the  diode.   Thus,  if  it  is  assumed  that  the  expression 

for  tunneling  current  remains  valid  in  the  presence  of 

impact  ionization  multiplication  (which,  in  light  of  the 

current  densities  encountered,  is  a  very  reasonable 

i 
assumption)  then  M   is  simply  the  ratio  of  the  measured 

value  of  reverse  current  to  the  extrapolated  value  of 

the  tunneling  current. 

Thus,  if 


S.,        2qI.M  (NR) 

l0g2qT^  =  l0g[ Iql ]  =  l0g[M  (NR)]        (3'8) 

i 
is  plotted  vs.  log  M  ,  as  done  in  Figures  29  to  32  for 

diodes  1N749  to  1N752,  a  straight  line  with  a  slope  of 

three  should  be  obtained.   It  is  clear  that,  with  the 

exception  of  the  1N752  diode  for  which  the  agreement 

between  experimental  results  and  theory  is  reasonable, 

considerable  discrepancies  exist.   Moreover,  the  discrep- 

S,d 
ancies  are  serious  because  the  measured  values  of  ,  J: 

2qlt 

are  larger  than  the  theoretical  one,  which  represents 
the  upper  limit  attainable  in  impact  ionization  multipli- 
cation for  which  electrons  and  holes  have  the  same  ioniza- 
tion coefficient.   In  addition,  at  low  values  of  multipli- 
ed 
cation  the  rate  of  increase  of  ~  ^   also  exceeds  the 

2qlt 


theoretical,  and,  therefore, the  upper  limit  value. 

At  this  point  it  might  be  argued  that  these  discrep- 
ancies arise  from  the  unequal  ionization  rates  of  electrons 
and  holes  in  silicon.   However,  if  the  latest  ionization 
rate  data  of  van  Overstraeten  and  de  Man  (1970) ,  shown  in 
Figure  33, are  extrapolated  to  the  field  strengths  encountered 
in  diodes  used  in  this  study,  the  extrapolated  electron 
ionization  rates  are  not  substantially  greater  than  the 
extrapolated  hole  ionization  rates. 

For  the  sake  of  completeness,  however,  the  noise 
data  was  also  analyzed  in  terms  of  Mclntyre's  (1966)  theory 
of  noise  current  spectral  density  for  impact  ionization 
multiplication  in  diodes  for  which  the  electron  and  hole 
ionization  rates  are  not  equal.   According  to  Mclntyre, 
if  the  carriers  initiating  the  multiplication  are  electrons, 
the  low  frequency  noise  current  spectral  density  is  given  by 


S.  =  2ql  (M')3[l-(l-k)  (^-^)2]  (3.9) 

r  M 


v/h 


ere  I   is  the  electron  current  initiating  the  multipli- 

g 
cation,  which  in  this  case  is  the  tunneling  current,  k=— 

is  the  ratio  of  the  hole  ionization  rate,  B,  to  the  electron 

■ 
ionization  rate,  a,  and  M   is  the  DC  current  multiplication 

,  I 
factor,  M  =-=^-    .   Since  3  and  a  depend  strongly  on  the 

electric  field,  this  equation  is  strictly  cor-pct  only  if 

an  appropriate  "average"  value  of  k  is  used  (Hclntyre ,  1973) 
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However,  for  the  purposes  of  this  discussion  this  is  not 

an  important  consideration. 

Si 

According  to  Mclntyre's  (1966)  theory,  ?     is  given  by 


i   .  /v,'  s3ri  „  ,_,  ,M  -1N2 


(M  )J[l-(l-k)  C-^r^)     ]  (3.10) 

Zqit  M 

which  is  plotted  for  several  values  of  k  in  Figure  34. 

Comparison  of  Figure  34  with  the  experimental  values  of 

S. 
1    of  Figures  29  to  31  indicates  that  a  closer  agreement 

between  theory  and  experiment  could  be  obtained  only  if 

the  hole  ionization  rate  were  substantially  larger  than 

the  electron  ionization  rate.   In  light  of  the  measured 

and  extrapolated  values  of  ionization  rates  of  Figure   33 

this,  however,  does  not  seem  very  likely. 

Moreover,  it  will  be  shown  in  Chapter  V  that  if  the 
measured  and  extrapolated  values  of  ionization  rates  of 
Figure  33  are  used  in  the  conventional  calculations  of 
DC  multiplication,  the  results  obtained  indicate  that 
multiplication  should  occur  for  any  non-zero  value  of 
reverse  bias.   This  implies  that  the  noise  ratio  should 
be  greater  than  unity  for  all  non-zero  values  of  reverse 
voltage,  which  is  in  direct  contradiction  with  the 
experimental  data  of  Figures  15  to  21. 

The  existing  theories  are,  therefore,  in  serious 
difficulty:   they  either  fail  to  correctly  predict  the 


s.  , 

magnitude  and  functional  dependence  of       on  M  ,  or 

2qlt 

require  an  unreasonable  ratio  of  ionization  rates  to 
possibly  achieve  an  acceptable  agreement  of  theory  and 
experiment;  they  fail  to  correctly  predict  the  magnitude 
of  the  experimentally  observed  DC  multiplication;  and 
they  completely  fail  to  predict  the  existence  of  a  reverse 
voltage  threshold  for  the  onset  of  impact  ionization 
multiplication.   Any  theory  of  reverse  breakdown  transition 
from  tunneling  to  impact  ionization  multiplication  must  be 
able  to  eliminate  all  of  these  discrepancies,  while 
remaining  consistent  with  all  other  observations  of  break- 
down phenomena  in  silicon  p-n  junctions. 

One  of  these  observations,  made  by  many  workers 
investigating  carrier  multiplication  in  semiconductors, 
is  the  existence  of  an  energy  threshold  for  pair  production 
(Shockley,  1961;  Moll  and  van  Overstraeten ,  19G3)  :   the 
agent  attempting  to  generate  a  carrier  pair  must  be  capable 
of  supplying  at  least  a  certain  minimum  amount  of  energy 
if  it  is  to  succeed  in  producing  a  hole-electron  pair. 
The  latest  values  of  threshold  energies  for  impact  ioniza- 
tion in  silicon,  as  determined  in  ionization  rate  measure- 
ments and  first  reported  by  Moll  and  van  Overstraeten  (1963) 
are  1.8±0.1  eV  for  electrons  and  2.4±0.1  eV  for  holes. 
Since  the  noise  data  of  Figures  15  to  21  indicate 
that  carrier  multiplication,  as  judged  from  the  departure 


of  the  noise  ratio  from  unity,  always  sets  in  at  a  well- 
defined  value  of  reverse  voltage  that  is  almost  independent 
of  the  diode  breakdown  voltage  (if  the  breakdown  voltage 
is  defined  as  that  value  of  reverse  bias  which  gives  rise 
to  some  constant  value  of  reverse  current  density) ,  the 
author  contends  that  the  onset  of  multiplication  is  governed 
by  the  magnitude  of  the  threshold  energy  for  ionization. 
To  substantiate  this  contention  quantitatively,  the 
following  multiplication  threshold  model,  described  in 
terms  of  Figure  35,  which  shows  the  electric  field  and 
energy  band  diagrams  for  a  reverse  biased  abrupt  p-n 
junction,  is  proposed.   Since  the  tunneling  probability 
is  strongly  dependent  on  the  electric  field,  most  of  the 
tunneling  current  originates  at  x=0.   At  x=x,  the  tunneling 
electrons  become  free  carriers,  absorbing  1.12±E   eV,  where 
±E   corresponds  to  the  emission  or  absorption  of  one  phonon 
whose  energy  may  be  17.9  meV  (TA) ,  43.7  meV  (LA),  53.2  meV 
(LO) ,  or  58.5  meV  (TO),  (Logan,  1969).   Continuing  towards 
x_  the  electrons  acquire  energy  from  the  electric  field, 
and  at  x„  those  which  have  made  the  fewest  number  of  colli- 
sions with  the  lattice  are  sufficiently  energetic  to 
participate  in  ionizing  collisions.   The  onset  of  impact 
ionization  multiplication,  therefore,  occurs  when  x.=x  . 
Denoting  the  reverse  voltage  at  which  this  occurs  by  V  ,  ,  , 
the  abrupt  junction  approximation  yields 


xn      (1.12±Er)+(1.8±0.1) 

x  +x   =        v — T+^ (3.11) 

n   p  rth  v 


where 


x  N.-N_ 

n  _   A   D 

x  +x  N 

n   p  A 


(3.12) 


Consequently 


rth 


[(1.12  +  E  )+(1.8±0.1)]  (_£-_) -<j>    .     (3.131 


Table  II  compares  the  reverse  threshold  voltaqes,  V   ,  , 

^     rth' 

computed  using  the  doping  densities  of  Table  I  and 
Er=0.06  eV,  with  the  reverse  threshold  voltages  obtained 
from  Figures  18  to  21.   In  light  of  the  approximations 
made   and  the  experimental  uncertainties  involved  in 
determining  the  doping  densities,  they  are  in  good  agreement 

The  holes  produced  in  ionizing  collisions  between 
x2  and  xn  travel  in  opposite  direction,  but  they,  too, 
ionize  only  after  they  have  acquired  at  least  the  necessary 
threshold  energy  for  ionization.   The  electrons  produced 

in  hole  ionizing  collisions  then  travel  towards  x  ,  and 

n 

after  acquiring  at  least  1.8±0.1  eV  participate  in  more 
ionizing  collisions.   This  process  repeats  itself,  but 
since  the  number  of  ionizing  collisions  decreases  with 
each  successive  set  of  ionizations,  the  process  converges 
to  a  finite  value  of  DC  multiplication. 
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TABLE  II 

REVERSE  VOLTAGES  AT  THE  ONSET  OF 

IMPACT  IONIZATION  MULTIPLICATION: 

COMPARISON  OF  DATA  WITH  THEORY 

DIODE  TYPE    EXPERIMENTAL  REVERSE    COMPUTED  REVERSE 
THRESHOLD  VOLTAGE      THRESHOLD  VOLTAGE 

1N749  2.65±0.1  V  2.99  -  3.45  V 

1N750  2.7+0.15  V  2.39  -  2.78  V 

1N751  2.7±0.15  V  2.31  -  2.71  V 

1N752  2.7±0.15  V  2.27  -  2.65  V 


In  light  of  this  model  it  is  not  difficult  to  see 
why  the  existing  theories  could  not  explain  the  experi- 
mental data:   none  of  the  theories  take  the  influence 
of  the  threshold  energies  for  ionization  explicitly  into 
account.   The  DC  multiplication,  which  will  be  discussed 
in  more  detail  in  Chapter  V,  has  conventionally  been 
calculated  by  integrating  the  ionization  rates  a(x)  and 
6(x)  over  the  junction  SCR  in  accordance  with  the  solution 
of  the  differential  equation 


dJ  (x)    dJ  (x) 


=  a(x)J  (x)+6(x)J  (x)        (3.14) 


dx       dx  n         p 

modeling  the  process  (Moll  and  van  Overstraeten,  1963; 
Lee   et  al.,  1964;  van  Overstraeten  and  de  Man,  1970). 

J  (x)  and  J  (x)  denote,  respectively,  the  electron  and 
n         p 

hole  current  densities  as  a  function  of  position  in  the 
junction  SCR.   This  equation  does  not  explicitly  consider 
the  effect  of  threshold  energies  on  the  multiplication 
process  because  it  claims  that  the  differential  increase 
in  current  densities  at  point  x  is  proportional  to  the 
current  densities  at  point  x,  which  does  not  take  into 
account  a  "distance  delay"  over  which  the  carriers  are 
incapable  of  ionizing  because  they  are  acquiring  the 
necessary  threshold  energy.   For  diodes  which  break  down 
at  several  hundred  volts  of  reverse  bias  this  "distance 
delay"  is  very  small  in  comparison  to  the  width  of  the 


junction  SCR,  and  the  approximation  involved  in  neglecting 
it  is  very  good.   For  diodes  breaking  down  at  several 
volts  of  reverse  bias,  on  the  other  hand,  the  "distance 
delay"  could  constitute  the  entire  width  of  the  diode  SCR. 
The  approximation  of  neglecting  it  in  such  cases  and 
employing  equation  3.14  is  totally  unacceptable  -  it  leads 
to  calculation  of  finite  DC  multiplication  when  no  multi- 
plication is  experimentally  observed. 

Application  of  conventional  calculations  of  DC  multi- 
plication to  very  low  breakdown  voltage  diodes   thus  leads 
to  the  conclusion  that  multiplication  must  be  taking  place 
in  these  junctions  because  the  electric  fields  in  them 
are  very  strong.   The  author,  on  the  other  hand,  contends 
that  if  the  electric  fields  in  these  junctions  are  very 
strong,  multiplication  will  occur,  but  only  if  the 
carrier  energies  exceed  the  threshold  energies  for  ioniza- 
tion.  This  is  precisely  why  no  impact  ionization  multi- 
plication occurs  in  tunnel  diodes,  in  spite  of  the 
extremely  strong  electric  fields  existing  in  these  junctions 

Because  Mclntyre's  theory  of  noise  in  avalanche 
diodes  employs  equation  3.14  in  its  derivation,  it  is 
now  clear  why  his  theory  is  not  capable  of  predicting 
the  multiplication  noise  in  diodes  employed  in  this  study. 
Moreover,  Tager's  theory  is  also  unacceptable  because 
it  does  not  incorporate  in  its  derivation  an  upper  limit 
on  the  number  of  ionizing  collisions  a  carrier  can  undergo 


in  its  transit  across  the  diode  SCR.   Such  a  limit  must 
be  imposed  if  threshold  energies  for  ionization  are 
taken  into  account. 

Not  explicitly  incorporating  the  influence  of  the 
threshold  energies  for  ionization  in  the  existing  theories 
of  noise  and  DC  multiplication,  therefore,  renders  them 
inapplicable  to  the  diodes  employed  in  this  study.   The 
next  two  chapters  will  attempt  to  resolve  these  difficulties, 


CHAPTER  IV 
DETAILED  ANALYSIS  OF  NOISE  DATA 

Figures  29  to  32  indicate  that  if  the  DC  multiplication 

i 
factor,  M  ,  is  defined  conventionally  by 

m'  =  ^  (4.D 

it  is  impossible  to  obtain  a  unique  slope  which  would  define 
the  exponent,  x,  in  Tager's  (1965)  expression  for  the  low 
frequency  noise  current  spectral  density 

S.  =  2qIt(M')X     .  (4.2) 

The  author  suspected  that  this  difficulty  arose  because 
the  DC  multiplication  factor  was  not  appropriately  defined. 
This  can  be  seen  most  clearly  when  M  =  1,  which  is  per- 
fectly acceptable  mathematically,  but  is  physically  meaning- 
less, for  a  multiplication  factor  of  unity  means  that  there 
is  no  carrier  multiplication.   Thus,  it  was  the  author's 
contention  that  the  DC  multiplication  factor  should  be 
defined  by 

I   =  It+MI.   ,  or  (4.3) 


M  =  -£  -  1  .  (4.4) 

t 
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This  definition  reduces  to  the  conventional  one  for  large 
values  of  M,  but  is  more  appealing  physically  for  it  allows 
M  to  go  to  zero  when  there  is  no  carrier  multiplication. 
However,  if  M  is  defined  by  M  =  (I  /I  )-l,  the  expression 
for  the  noise  current  spectral  density  has  to  be  reformulated, 
for  consistency,  to 

Si  =  2qlt[l+f (M)]  (4.5) 

where  f (M)  goes  to  zero  when  M  goes  to  zero. 

The  results  of  plotting  f (M)  vs.  M  for  diodes  1N749 

to  1N752  are  shown  in  Figures  36  to  39.   These  graphs 

contain  several  interesting  features.   First,  even  though 

all  data  points  do  not  lie  on  one  straight  line,  it  is 

possible  to  fit  most  of  them  along  two  or  three  straight 

lines,  the  slopes  of  which  increase  monotonically  with 

increasing  reverse  voltages.   Second,  the  departures  of 

data  points  from  lines  of  lower  slopes  to  lines  of  higher 

slopes  occur   for  all  diodes  at  roughly  the  same  (although 

somewhat  increasing)  value  of  reverse  voltage.   The  values 

of  reverse  voltages  at  which  these  transitions  set  in,  as 

defined,  whenever  possible,  by  the  intersection  of  the  two 

slopes,  are: 

Diode  type  Onset  of  1st  Onset  of  2nd 

transition  transition 

1N749  3.95  V  4.35  V 

1N750  4.1  V  4.65  V 

1N751  4.4V  4. 95  V 

1N752  4.7  V  5.1  V 


Finally,  even  though  the  slope  of  any  line  tangent  to  the 
data  points  does  not  exceed  its  maximum  allowed  value  of 
three,  as  did  the  tangents  to  the  data  points  of  Figures 
36  to  39,  it  is  now  possible  to  draw  tangent  lines  whose 
slopes  approach  unity  --  which  appears  to  contradict  the 
established  results  that,  according  to  the  old  definitions, 
the  values  of  the  slopes  should  never  be  less  than  two. 

All  of  these  observations  can  be  understood  in  terms 
of  the  conduction  model  incorporating  the  carrier  threshold 
energy  for  ionization,  as  proposed  in  the  preceding  chapter. 
Since  the  amount  of  noise  produced  in  impact  ionization 
multiplication  increases  as  the  maximum  number  of  ionizing 
collisions  per  carrier  transit  across  the  SCR  increases, 
and  eventually  approaches  the  upper  limit  value  given  by 
Tager's  (1965)  expression  for  the  case  of  unlimited  number 
of  ionizing  collisions  on  a  carrier  transit  across  the  SCR, 
the  breaks  in  the  straight  lines  in  Figures  3  6  to  3  9  and 
the  monotonic  increase  in  slopes  are  indicative  of  transi- 
tions from  a  lower  to  a  higher  maximum  number  of  ionizing 
collisions  that  a  carrier  can  undergo  in  a  transit  across 
the  diode  SCR.   In  particular,  if  the  built-in  voltage  of 
the  junctions  is  taken  into  account,  and  an  electron  threshold 
energy  of  1.8±0.1  eV  is  assumed,  as  given  by  Moll  and  van 
Overstraeten  (1963) ,  then  it  becomes  clear  that  the  onset 
of  the  first  break  in  the  slopes  of  Figures  36  to  39  is 
due  to  a  transition  from  a  maximum  of  two  to  a  maximum  of 


three  ionizing  collisions  per  electron  transit  across 

the  diode  SCR.   This  is  easily  verified  by  dividing 

1.8±0.1  eV  into  the  total  diode  potential,  the  sum  of 

the  built-in  potential  listed  in  Table  1,  and  the  potential 

at  the  onset  of  the  first  slope  transition,  listed  above. 

The  results  are: 

Diode  type       Total  potential  divided  by  1.8±0..  1  eV 

1N749  2.92  -  2.62 

1N750  3.01  -  2.69 

1N751  3.19  -  2.85 

1N752  3.36  -  3.01 

When  the  procedure  is  reversed  and  the  total  potential 
at  the  onset  of  the  first  slope  transition  is  divided 
by  three,  the  following  electron  threshold  energies  for 
ionization  are  obtained: 

Diode  type       Electron  threshold  energy  for  ionization 

1N749  1.66  eV 

ln750  1.7  eV 

1N751  1.8  eV 

1N7  52  1.9  eV 

This  is  in  excellent  agreement  with  Moll  and  van  Over- 
straeten's  (1963)  value  of  1.8±0.1  eV.   The  apparent  increase 
in  the  threshold  energy  for  ionization  with  increasing 
diode  breakdown  voltage  is  also  reasonable.   Both  Moll  and 
van  Overstraeten' s  (1963)  and  these  values  include  the 
"intrinsic"  threshold  energy  as  well  as  energy  lost  in 


phonon  collisions.   Since  the  rate  of  energy  acquisition  with 
distance  is  smaller  in  higher  voltage  breakdown  diodes  (the 
electric  fields  are  weaker) ,  but  the  rate  of  energy  losses 
is  the  same  (because  it  depends  on  mean  free  path  between 
phonon  collisions,  which,  presumably,  is  not  strongly 
dependent  on  field  strength) ,  it  can  be  expected  that  an 
increasingly  larger  number  of  phonon  collisions  would  occur 
between  successive  ionizations,  resulting  in  increasingly 
higher  observed  threshold  energies  in  higher  breakdown 
voltage  diodes. 

An  analogous  analysis  can  be  made  for  the  second  slope 
transition  in  Figures  36  to  39.   Since  the  threshold  energy 
for  ionization  for  holes  is,  according  to  Moll  and  van 
Overstraeten  (1963) ,  greater  than  the  threshold  energy  for 
ionization  for  electrons,  the  second  slope  transition  may 
be  suspected  to  arise  from  a  transition  from  a  maximum  of 
two  to  a  maximum  of  three  ionizing  collisions  per  hole 
transit  across  the  diode  SCR.   When  the  total  potential 
at  the  onset  of  the  second  slope  transition  is  divided 
by  three,  the  following  values  of  threshold  energy  for 
ionization  for  holes  are  obtained: 

Diode  type       Hole  threshold  energy  for  ionization 

1N749  1.79  eV 

1N750  1.89  eV 

1N751  1.99  eV 

1N752  2.04  eV 


3  3 


These  values  are  considerably  smaller  then  the  2.4+0.1  eV 

obtained  by  Moll  and  van  Overstraeten  (1963) .   Thus  it  may 

be  suspected  that  the  second  slope  transition  is  due  to  a 

transition  from  a  maximum  of  one  to  a  maximum  of  two  ionizing 

collisions  per  hole  transit  across  the  SCR.   This,  however, 

results  in  hole  threshold  energies  that  are  substantially 

higher  than  those  given  by  Moll  and  van  Overstraeten  (1963) : 

Diode  type        Hole  threshold  energy  for  ionization 

(assuming  transition  from  a  maximum  of 
one  to  a  maximum  of  two  ionizations  per 
transit) 

1N749  2.69  eV 

1N750  2.84  eV 

1N751  2.99  eV 

1N752  3.06  eV 

Moreover,  if  the  hole  threshold  energies  obtained  under  this 
assumption  are  used  in  the  simulation  of  DC  multiplication, 
a  serious  difficulty  arises  which  does  not  occur  when  the 
hole  threshold  energy  of  1.9±0.15  eV  is  used.   As  a  con- 
sequence, it  is  suggested  that  the  value  of  2.4±0.1  eV  is 
probably  too  high. 

An  immediate  application  of  this  interpretation  of 
noise  data  can  be  made  in  conjunction  with  the  conduction 

model  proposed  in  the  preceding  chapter  and  the  — 5-  vs.  V 

C^ 

data  to  compute  the  doping  densities  on  the  n  and  p  sides 
of  the  junctions.   With  the  help  of  Figure  35,  the  analysis 
proceeds  as  follows.   At  the  onset  of  impact  ionization  multi- 
plication 


V  +V     N  -N 

where  V  =  E  /q,  V   =  E  /q  and  E   is  the  bandgap  energy 

and  E   is  the  electron  threshold  energy  for  ionization. 

Moreover,  from  — j   vs.  V  data 
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data,  and  V  ^,  and  V  can  be  determined  from  the  noise  data, 
rth      e 

N   and  N   can  be  calculated.   A  comparison  of  N^  and  N_ 
DA  v  D       A 

obtained  in  this  manner  with  the  values  of  N   and  N  given 
in  Table  I  is  shown  in  Table  III.   The  agreement  is  quite 
good,  which  lends  further  credence  to  the  proposed  conduction 
model  and  the  subsequent  interpretation  of  noise  data. 
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Moreover,  if  the  proposed  conduction  model  is  truly 
valid,  then  a  suitable  derivation  incorporating  the  essential 
features  of  the  model  should  yield  expressions  for  the  mag- 
nitude of  the  noise  current  spectral  density  which  should  be 
in  good  agreement  with  experimentally  observed  results. 
From  the  point  of  view  of  noise  as  a  statistical  fluctuation 
phenomenon,  the  most  important  feature  of  the  proposed  con- 
duction model  is  the  limit  it  imposes  on  the  maximum  number 
of  ionizing  collisions  a  carrier  can  experience  in  one 
transit  across  the  diode  SCR.   Consequently,  any  derivation 
of  the  noise  current  spectral  density  must  incorporate  this 
constraint.   This  complicates  the  problem  considerably  for 
it  requires  that  the  multiplication  process  be  examined 
on  basis  of  individual,  successive  transits  of  carriers 
across  the  diode  SCR,  rather  than  on  monitoring  the  end 
product  of  the  multiplication  process  at  the  device  terminals 
as  was  done  by  Tager  (1965)  and  Mclntyre  (1966) . 

To  obtain  some  insight  into  how  the  analysis  of  multiple 
ionizing  collisions  might  be  carried  out,  the  author  began 
by  considering  the  limiting  case  of  a  maximum  of  one  ionizing 
collision  per  carrier  transit  across  the  diode  SCR.   This 
is  more  than  just  an  exercise.   According  to  the  author's 
interpretation  of  Figures  36  to  38,  about  half  of  the 
data  points  correspond  to  multiplication  resulting  from 
no  more  than  two  ionizing  collisions  per  carrier  transit 
across  the  SCR.   Moreover,  even  when  up  to  two  ionizing 
collisions  are  energetically  possible,  the  probability  of 


the  second  collision  will  be  considerably  smaller  then 
the  probability  of  the  first  collision  because  the  space 
charge  regions  of  the  diodes  examined  in  this  study  are 
very  narrow,  whereas  the  mean  free  path  between  ionizing 
collisions,  as  quoted  by  previous  workers  (Moll  and  van 
Overstraeten,  1963;  Lee   et  al.,  1964)  is  relatively 
large.   Thus  the  approximation  of  no  more  than  one  ionizing 
collision  per  carrier  transit  across  the  diode  SCR  should 
yield  reasonable  results,  especially  for  the  case  of  the 
diode  with  the  narrowest  SCR,  the  1N74  9. 

The  analysis  of  the  noise  current  spectral  density 
is  best  carried  out  by  referring  to  Figure  35.   The  tunneling 
electrons,  originating  at  x=0  at  the  rate  of  n   per  second, 
are  accelerated  from  x=x,  to  x=x_,  and  undergo  ionizing 
collisions  from  x=x2  to  x=x  ,  producing  n,  electron-hole 
pairs.   The  n   tunneling  electrons  and  the  n.  ionization 
produced  electrons  are  swept  to  the  right  and  out  of  the 
SCR,   The  n1   holes  travel  to  the  left  and,  after 
acquiring  the  necessary  threshold  energy,  undergo  ionizing 
collisions,  producing  n   electron-hole  pairs.   The  (n  +n  ) 
holes  are  then  swept  to  the  left  and  out  of  the  SCR,  while 
the  n2  electrons  travel  to  the  right  and,  after  acquiring 
the  necessary  threshold  energy,  ionize,  producing  n^  electron- 
hole  pairs.   This  process  continues,  so  that  the  total 
number  of  carriers,  N,  collected  at  the  device  terminals 
becomes 


I  n      •  (4.11) 

i=0  x 


Letting  n-^a^  ,  n  =b_n1=b2a1n  ,  n^a^^a^^.^.,  etc., 
where  the   denotes  averages  and  a^  and  b^  represent  the 
probability  of  ionization  for  electrons  and  holes,  respectively, 
on  the  i-th  transit  across  the  diode  SCR, 


N  =  nn      E    n  a.   _,b.       ;  a  =1,  b  =1  .    (4.12) 
0  i=0  j=0  D  odd  3  even     0      0 

This  converges  to  a  finite  number  because  a.<l  and  b.<l 
since  no  more  than  one  ionizing  collision  per  carrier  transit 
across  the  SCR  is  allowed.   A  reasonable  simplification 
which  can  be  made  at  this  point  is  that  a . =a  and  b.=b.   This 
means  that  the  probabilty  of  ionization  differs  for  holes  and 
electrons,  but  remains  the  same  on  all  transits.   In  this  case 


N  =  n0(l+a+ab+a2b+a2b2+  ...) 


=  nn(^S-)  (4.13) 

u  1-ab 

and  the  DC  multiplication  factor,  M  =  —  -  1,  becomes 
m   i(l+b) 

1-5B        *  (4-14) 

The  low  frequency  spectrum  of  the  current  fluctuations 
resulting  from  the  multiplication  process  can  be  obtained  from 


SM  =  2var  N  (4.15) 

N 


where  SN  denotes  the  magnitude  of  the  spectrum  of  N,  the 

number  of  events  occurring  during  a  unit  time  interval 

(van  der  Ziel,  1970).   Since  the  diode  current,  I  ,  is 

r 

given  by  I  =qN, 

2 

Sj   =  2q  var  N  (4.16) 

r 

where 


~2   -2 

var  N  =  N   -N 


=  (  Z  n.)2  -  (  I    n.) 
i=0  1      i=0  X 


=  T.        I    (n.n.-n.n.)      .  (4.17) 

i=0  j-0   1  ^   i  1 

In  order  to  evaluate  S    it,  therefore,  remains  to  evaluate 

r 

n.n  .   This  is  accomplished  by  representing  n.  by 


0  b .  ;  i  even 

ni  =  rf1alrb2ra3rb4r  •"   "  (4'18) 

a .  ;  i  odd 
lr 

where  a.   and  bmr  are  either  1  or  0  (depending  on  whether  or 

not  an  electron  or  hole  ionizes  on  the  j-th  or  m-th  transit 

across  the  SCR)  and  ensemble  averaging  over  the  products 

n. n. .   This  analysis  is  discussed  in  detail  in  Appendix  A. 

The  resulting  expression  for  var  N  is  given  by 

var  N  =  n  {l±3i±3ib+£2l}  {419) 

u   (1-55)^ 


and,  therefore, 


_  2-  fl+3a+3ab+a2bi 

SI   =  2q  n0f =~2 

r        °   (1-abr 


=  2gI  {l+35+3ab+a2b|        .  (4>2o; 

t     (1-abp 

If,  for  the  moment,  it  is  assumed  that  the  probability 
of  ionization  is  the  same  for  holes  as  it  is  for  electrons, 
a=b,  then 


(4.21) 


1-a 


and 


ST   =  2ql.  a+a\  .  (4.22) 

Jr       t(l-a)/ 

Eliminating  a, 

SI   =  2qIt(1+3M+2M2)  (4.23) 

r 

which  is,  indeed,  in  the  previously  postulated  form 

Sj   =  2qIt[l+f(M)]     .  (4.24) 

This  result  provides  a  more  formal  justification  of  the 
previously  proposed  redefinitions  of  the  DC  multiplication 
factor  and  the  functional  form  of  S.,  and  lends  greater 
credence  to  the  interpretation   of  noise  data  based  on  the 
new  definitions. 


41 


When  equation  4.20   is  rewritten  in  the  form 

S    =  2ql  [l+f(5,b)]  (4.25) 

r 

and  f  (a,b)  is  plotted  as  a  function  of  M  for  various  values 

of  k=— ,  Figure  4  0  is  obtained,  which  contains  the  somewhat 

a 
surprising  result  that  f (a,b)  is  rather  weakly  dependent  on  k. 

This  result  is  exploited  in  the  derivation  of  the  magnitude 

of  the  noise  current  spectral  density  for  the  case  of  up 

to  two  ionizations  per  carrier  transit  across  the  diode  SCR. 

Moreover,  since 

Sj   =  2qlt[l+f (a,b)]  *  2ql  (1+3M+2M2)  (4.26) 

over  a  fairly  wide  range  of  k=^  for  low  values  of  M,  it 

a 
demonstrates  that  the  slope  of  f (M)  vs.  M  can,  indeed, 

approach  unity,  as  it  does  for  the  data  of  Figure  36. 

To  compare  the  results  of  this  analysis  with  experiment, 

2 
Figure  41  shows,  along  with  f (M)=3M+2M  plotted  versus  M, 

the  experimental  data  for  the  1N74  9  diode.   Although  the 

theoretical  curve  shows  the  same  functional  dependence  on 

M  as   do   the  experimental  data,  it  is  about  a  factor  of 

three  too  low.   However,  in  view  of  the  fact  that  the 

derivation  of  the  theory  is  strictly  statistical,  taking 

no  account  of  any  junction  parameters,  this  result  was 

most  encouraging. 

At  this  point  the  author  suspected  that  the  discrepancy 

between  theory  and  data  resulted  from  dismissing  the  possibility 


of  the  second  and  third  ionizations,  and  proceeded  to 
extend  the  theory  to  more  than  one  ionization  per  carrier 
transit  across  the  SCR.   Although  this  is  easily  done  in 
principle,  it  becomes  extremely  tedious  algebraically 
because  the  analysis  must  take  into  account  the  various 
branching  possibilities,  and  the  correlation  between  the 
branches.   In  fact,  the  extension  of  the  theory  to  the 
case  of  up  to  two  ionizations  per  carrier  transit  across 
the  SCR  is  a  composite  of  an  infinite  number  of  the  up-to- 
one-ionization  analyses.  In  order  to  restrict  the  algebraic 
details  to  manageable  size,  the  author  chose  to  analyze 
the  case  of  up  to  two  ionizations  per  carrier  transit  under 
the  additional  restriction  that  the  probability  of  the 
second  ionization  is  much  smaller  than  the  probability  of 

the  first  ionization,  and,  therefore,  that  terms  containing 

-2 

a»,  where  a2  is  the  probability  of  the  second  ionization, 

could  be  neglected.   It  was  felt  that  this  would  be  a 
reasonable  approximation  for  the  noise  data  of  Figures  3  6 
to  39  below  the  first  slope  transition. 

The  analysis  of  this  case  under  the  additional  approxi- 
mation that  the  electron  and  hole  ionization  probabilities 
are  the  same  (as  justified  by  DC  calculations  of  Chapter  V 
and  by  the  relative  insensitivity  of  the  results  of  Figure  40 
to  k=— )  is  given  in  Appendix  B.   The  algebraically  tedious 

1  *2 

results,  evaluated  for  ^—  =  0.3,  which  is  the  highest  ratio 


of  the  probability  of  the  second  ionization  to  the  probability 

of  the  first  ionization  that  the  author  felt  could  be  assumed 

under  the  approximation  that  a2<<ai'  are  compared  in  Figure  4  2 

with  the  results  of  the  up  to  one  ionization  per  transit 

theory,  evaluated  for  k=-^=l.   Although  the  two  curves  inter- 

a 
sect  giving  rise  to  a  slope  transition,  in  accordance  with 

the  previous  interpretation  of  the  slope  transitions  in 
the  experimental  data,  the  results  are  somewhat  disappointing 
in  that  f  (a,,a2)  is  not  sufficiently  larger  to  give  sub- 
stantially better  agreement  between  theory  and  the  noise 
data  of  Figures  36  to  39. 

A  careful  examination  of  the  various  assumptions  and 
approximations  made  in  the  theoretical  analysis  revealed 
that  the  source  of  the  final  discrepancy  between  the  theoret- 
ical results  and  the  experimental  data  lay  in  the  simplifying 
approximation  that  the  carrier  ionization  probability  is  the 
same  on  all  transits.   A  typical  plot  of  the  number  of 
carriers  generated  on  a  given  transit  per  one  tunneling 
electron,  obtained  from  DC  calculations  for  the  1N752  diode 
operating  at  reverse  bias  of  5.2  volts,  and  shown  in 
Figure  43,  demonstrates  that  the  ionization  probability 
for  electrons  on  the  first  transit  across  the  diode  SCR 
is  considerably  smaller  than  the  ionization  probability 
for  carriers  traversing  the  SCR  on  succeeding  transits. 
This  may  be  deduced  by  comparing  the  ratio  of  consecutive 


terms  for  transit  numbers  greater  than  about  five  which, 
in  the  case  of  Figure  42  is  0.88,  with  the  value  of  the 
first  term, 0.516.   If  the  ionization  probability  were 
identical  for  all  transits,  the  value  of  the  first  term 
would  have  been  0.88.   In  retrospect,  it  is  easy  to  see 
why  the  ionization  probability  on  the  first  transit  is 
less  than  on  succeeding  transits.   The  electrons  under- 
going ionizations  on  the  first  transit  come  from  tunneling, 
which,  in  terms  of  Figure  35,  means  that  they  do  not  become 
free  carriers  until  x=x, .  Consequently,  they  ionize  over 
a  much  smaller  portion  of  the  SCR  than  do  the  succeeding 
holes  or  electrons  originating  as  free  carriers  closer 
to  the  edges  of  the  SCR. 

To  determine  what  effect  the  smaller  ionization 
probability  on  the  first  transit  would  have  on  the  noise 

theory,  the  derivation  of  the  up-to-one-ionization-per- 

-  -   *1 
SCR-transit  theory  was  repeated  for  the  case  of  a=b,   ^-=r<l, 

a 

where  a.  is  the  ionization  probability  on  the  first 

transit.   The  result,  derived  in  detail  in  Appendix  A, 

is  given  by 

S.  =  2qIt(l+3M+2M2/r)  (4.27) 

clearly  indicating  that  theoretically  predicted  noise 
output  will  increase  for  r<l. 

This  analysis  was  also  repeated  for  the  up  to  two 
ionizations  per  SCR  transit  theory.   Indeed,  this  is  the 
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derivation  that  appears  in  Appendix  B.   Because  the  results 
are  algebraically  tedious,  they  are  evaluated  numerically 
once  the  values  of  a^,  a^,  and  a"2  have  been  determined. 
Here,  ai;L  denotes  the  probability  of  the  first  ionization 
on  the  first  transit,  and  a^  and  a„  are  the  first  and 
second  ionization  probabilities  on  succeeding  transits. 
For  the  range  of  voltages  where  the  up  to  two  ionizations 
theory  is  valid,  the  probability  of  second  ionization  on 
the  first  transit  is  zero.   The  details  of  evaluating 

all'  al  and  a2  from  the  DC  multiplication  program  are 
discussed  in  Appendix  C. 

The  final  results  of  the  up  to  two  ionizations  per 
SCR  transit  theory  for  the  case  of  a11<a,  are  compared 
with  the  experimental  data  for  the  1N74  9  to  1N751  diodes 
in  Figures  44  to  46.   In  view  of  the  many  approximations 
made  in  the  derivation  of  the  theory,  the  agreement  is 
quite  good,  even  for  those  values  of  M  where  the  probability 
of  a  third  ionization  is  no  longer  zero.   In  the  case  of 
the  1N751  diode,  where  at  low  values  of  M  the  theory  predicts 
a  larger  noise  output  than  appears  to  have  been  observed, 
the  discrepancy  is  probably  due  to  incorrect  values  of  M, 
which  at  low  values,  can  be  subject  to  substantial  error. 

The  analysis  of  the  noise  data  was  terminated  at  this 
point  for  it  was  felt  that  even  though  an  up-to-three- 
ionizations-per-SCR-transit  theory  would  undoubtedly 


yield  better  agreement  between  theory  and  experiment,  it 
would  probably  offer  no  new  insights  into  the  details  of 
the  multiplication  mechanism,  and,  therefore,  in  view  of 
the  algebraic  tedium  required  for  its  derivation,  would 
not  be  worthwhile. 


CHAPTER  V 
CALCULATION  OF  DC  MULTIPLICATION 

As  has  already  been  mentioned  in  Chapter  III,  the 
calculation  of  DC  multiplication  in  avalanche  diodes  has 
conventionally  been  done  by  integrating  the  electron  and 
hole  ionization  rates,  a  and  6,  over  the  junction  SCR, 
in  accordance  with  the  solution  of  the  differential  equation 


dJ  (x)    dJ  (x] 
2 =  S 


a(x)J  (x)+B(x)J  (x)       (5.1) 


dx        dx  n         p 

modeling  the  process.   Since  this  equation  incorporates 
the  effect  of  the  threshold  energies  for  ionization 
implicitly  in  the  ionization  rates  and,  therefore,  is 
not  capable  of  modeling  those  segments  of  the  carrier 
trajectories  over  which  multiplication  does  not  occur 
because  the  carriers  do  not  have  sufficient  energy  to 
ionize,  it  is  not  applicable  to  calculation  of  DC 
multiplication  in  low  breakdown  voltage  diodes  in  which 
these  "inactive"  segments  of  carrier  trajectories  con- 
stitute a  very  sizeable  portion  of  the  diode  SCR. 

These  "inactive"  segments  could  be  taken  into  account 
in  the  above  equation  by  incorporating  "distance  delays," 
L  (x)  and  L  (x)  ,  so  that 


47 


(5.2) 
dJ  (x)    dj  (x) 
^—  =    gx    =  a(x)Jn[x-Ll(x)]+6(x)Jp[x+L2(x)] 

However,  if  L,(x)  and  L  (x)  are  taken  so  that 

x 

V   -  -q/   E(x)  dx  (5.3) 

e       x-L1(x) 

and 


x+L  (x) 
■qf         *    E(x)  dx  (5.4) 


where  V  and  V,  represent,  respectively,  the  threshold 

potentials  for  ionization  of  electrons  and  holes,  then 

the  magnitudes  of  J   and  J  must  be  increased  from  their 
n      p 

actual  values  at  x±L(x)  since,  in  general,  there  will  be 
some  carriers  with  non-zero  energies  which  will  not  need 
to  travel  the  full  distance  L(x)  to  acquire  the  necessary 
threshold  energy  for  ionization.   On  the  other  hand,  if  J 
and  J   are  left  unaltered,  then  L,  (x)  and  L?  (x)  need  to  be 
evaluated  as  suitable  averages  over  all  carrier  energies. 
Since  neither  correction  is  easy  to  make  at  the  outset,  and 
since  the  modified  differential  equation  would  still  need 
to  be  solved  numerically,  the  author  chose  to  abandon  this 
approach  in  favor  of  a  numerical  simulation  of  the  mul- 
tiplication process. 

The  simulation,  whose  program  listing  is  given  in 
Appendix  C  is  accomplished  as  follows.   An  abrupt  junction 


doping  profile  is  assumed,  and  the  electric  field  triangle 
is  cut  into  a  large  number  of  variable  width  segments  such 
that  the  integral  of  the  electric  field  is  the  same  over 
each  segment.   (In  the  results  shown  in  Figures  47  to  50, 

x+Ax 
dV  =  -q/   E(x)  dx  =  0.01V,  so  that  the  electric  field 
x 

triangle  was  cut  into  several  hundred  segments.)   The 

integral  of  the  ionization  rates,  which  represents  the 

probability  of  an  ionization,  and  the  fraction  of  the 

tunneling  current  originating  in  a  given  interval  are 

then  evaluated  approximately  for  each  interval  by 

VAxk  , 

dPn(k)  =  /   a(x)  dx  *  ^[a  (xk+Ax  )  +a  (xk)]  Axk  (5>5) 

VAxk 
dPh(k)  =  /    6(x)  dx  «  ^[6(xk+Axk)+6(xk)]Axk  (5<6) 

Xk 
and 

Xk+Axk   s' 

/      "FTxT  , 
x,    e      dx 


dlt(k)  - 


w     s 
/   e_F^"dx 


0 


l[e  E(xk+Axk)  +  e  E(xk)]AXk 

~~~? — " — ZZEZ — 

iz[e"fc(xk+Axk)  +  e_FTxy"]Axk 


(5.7) 


where  dP  (k)  and  dP,(k)  represent,  respectively,  the 

probability  of  ionization  for  electrons  and  holes  in 

the  k-th  interval,  and  dl  (k)  represents  the  normalized 

tunneling  current  originating  in  the  k-th  interval.   The 

multiplication*  process  is  then  simulated  by  following  the 

motion  of  electrons  and  holes  back  and  forth  across  the 

diode  SCR,  and  allowing  them  to  participate  in  ionizations 

after  they  acquire  their  respective  threshold  energies, 

which  can  be  kept  track  of  very  easily  by  counting  off 

the  appropriate  number  of  intervals.   Thus,  the  tunneling 

current  originating  in  the  k-th  interval  becomes  a  free 

V 
electron  current  in  the  (k+m)-th  interval,  where  m=-r~L 

which  may  then  ionize  starting  in  the  (k+m+n,)-th  interval, 

V  L 

where  n,=g^,  and  where  V  and  V  denote,  respectively,  the 

bandgap  potential  and  the  electron  threshold  potential  for 

ionization.   That  fraction  of  the  tunneling  current  which 

ionizes  in  a  particular  interval,  and  the  ionization  produced 

electrons  originating  in  that  interval,  are  then  not  allowed 

to  ionize  again  until  they  have  traveled  at  least  another 

n,  intervals.   The  tunneling  electrons  are  traced  in  this 

manner  across  the  entire  SCR,  while  the  holes  produced  in 

ionizations  occurring  in  a  given  interval  are  traced  in  the 

opposite  direction,  not  being  allowed  to  ionize  until  they 

vh 

have  traveled  at  least  n  =-3^  intervals,  where  V,  is  the  hole 
2.    dv  h 

threshold  potential  for  ionization.   The  fraction  of  the 

hole  current  ionizing  in  a  particular  interval,  along  with  the 


ionization  produced  holes  originating  in  that  interval, 
are  then  not  allowed  to  ionize  again  until  they  have 
traversed  at  least  another  n   intervals.   In  this  manner 
the  hole  current  is  traced  across  the  entire  SCR. 

This  procedure  is  alternately  repeated  for  electrons 
and  holes,  and  the  contribution  to  the  terminal  current 
due  to  ionizations  on  each  carrier  transit  across  the  SCR 
is  recorded.   When  the  contribution  on  the  p-th  transit 
is  less  than  0.005  of  the  contribution  on  the  first  transit, 
the  simulation  is  stopped  and  the  contributions  to  the 
terminal  current  from  all  transits  are  added.   The  same 
procedure  is  followed  for  the  multiplication  process  initiated 
by  the  holes  left  in  the  valence  band  when  the  electrons  tunnel  to 
the  conduction  band.   The  total  multiplication  is  then  the 
sum  of  the  contributions  from  the  tunneling  electron 
initiated  multiplication  and  the  typically  much  smaller 
hole  initiated  multiplication. 

The  simulation  thus  circumvents  both  of  the  previously 
discussed  difficulties  which  would  be  encountered  in  the 
differential  equation  approach  to  the  computation  of  DC 
multiplication.   It  incorporates  unambiguously  the  threshold 
energies  for  ionization  and,  by  keeping  track  of  the  past 
history  of  the  carriers,  automatically  computes  the  fraction 
of  the  total  current  crossing  a  given  interval  which  is 
sufficiently  energetic  to  ionize  in  that  interval.   This  is 
equivalent  to  automatically  adjusting  the  magnitudes  of 


J  [x-L  (x) ]  and  J  [x+L  (x) ]  to  account  for  the  fact  that 
not  all  carriers  need  to  travel  the  full  distance  L(x)  in 
order  to  acquire  the  necessary  threshold  energy  for  ionization. 

Unfortunately,  the  simulation  creates  a  different 
problem  which,  too,  is  difficult  to  resolve  at  the  outset. 
Namely,  it  assumes  that  when  an  ionizing  collision  occurs 
the  primary  carrier  and  the  secondary  carriers  are  scattered 
in  a  way  which,  on  the  average,  does  not  impart  any  excess 
energy  to  either  the  primary  or  the  secondary  carriers,  so 
that  all  carriers  will  be  accelerated  again  essentially 
from  rest.   Since  the  scattering  produced  by  an  ionizing 
collision  is  probably  not  isotropic,  this  assumption  may 
not  seem  very  reasonable.   However,  in  the  actual  calcula- 
tions the  non-isotropic  effects  should,  at  least  partially, 
cancel,  for  if  the  scattering  favors  the  direction  of  the 
ionizing  carrier,  then  the  carriers  of  the  same  polarity 
as  the  ionizing  carriers  will  not  need  to  accelerate  the 
full  distance  L(x)  to  acquire  the  threshold  energy,  whereas 
the  carriers  of  opposite  polarity  will  need  to  accelerate 
farther  than  the  distance  L(x)  to  acquire  their  threshold 
energy.   Since  the  simulation  is  performed  by  alternately 
tracing  the  motion  of  electrons  and  holes  across  the  SCR, 
these  effects  should  alternately  cancel  to  some  extent. 
However,  without  knowing  the  exact  details  of  the  ionization 
and  scattering  mechanism,  it  is  difficult  to  estimate  how 
much  cancelling  actually  takes  place. 
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Although  the  simulation  program  was  written  especially 
for  computation  of  DC  multiplication  occurring  in  low  break- 
down voltage  diodes,  where  it  is  necessary  to  take  the 
threshold  energies  for  ionization  explicitly  into  account, 
it  may  also  be  used  to  perform  DC  multiplication  calculations 
in  the  conventional  manner  simply  by  letting  the  threshold 
energies  for  ionization,  and  the  bandgap  energy,  approach 
zero.   In  actual  calculations  the  bandgap  and  threshold 
potentials  were  set  at  0.01  V.   The  ionization  rates 
employed  were  those  of  van  Overstraeten  and  de  Man  (1970) . 
From  the  results  obtained  for  diodes  1N750  to  1N752,  shown 
in  Figures  47  to  49  under  the  label  of  "conventional 
calculations,"  it  is  apparent  that  the  conventional  method 
of  performing  DC  multiplication  calculations  is  inadequate 
for  it  predicts  much  greater  multiplication  than  is  experi- 
mentally observed,  even  to  the  point  of  predicting  finite 
multiplication  when  no  multiplication  actually  occurs,  as 
evidenced  by  the  noise  data  of  Figures  15  to  21. 

Although,  as  shown  in  Figures  47  to  4  9  under  the 
label  of  "conventional  calculations,  including  tunneling," 
the  agreement  between  the  experimental  and  the  calculated 
results  improves  somewhat  when,  on  the  first  transit,  the 
bandgap  energy  for  tunneling  is  taken  into  account,  the 
calculated  values  of  DC  multiplication  are  still  considerably 
larger  than  what  is  experimentally  observed,  and  finite 
values  of  multiplication  are  still  calculated  when  no 
multiplication  actually  occurs. 


When  the  threshold  energies  for  ionization  are 
finally  taken  into  account,  the  calculated  values  of 
DC  multiplication  do  go  to  zero  when  no  multiplication 
is  observed  in  the  noise  data,  but  the  predicted  magnitudes 
of  the  DC  multiplication  are  smaller  than  those  experimentally 
observed.   This,  however,  is  understandable  for  the  integral 
of  the  ionization  rates  over  a  given  distance  represents, 
by  definition,  the  probability  of  ionization  over  that 
distance.   Therefore,  eliminating  large  portions  of  the 
diode  SCR  from  integration  by  explicitly  taking  the  threshold 
energies  for  ionization  into  account,  but  employing  conven- 
tional ionization  rates  derived  from  integration  over  the 
entire  SCR,  reduces  the  calculated  ionization  probabilities 
and,  therefore,  the  calculated  DC  multiplication  factor. 
It,  therefore,  seems  appropriate  to  increase  the  ionization 
rates  to  compensate  for  this  shortening  of  the  ionization 
rates'  integration  region.   Unfortunately,  this  is  not  a 
trivial  matter  since  the  functional  dependence  of  the 
ionization  rates  on  the  electric  field  implicitly  reflects, 
along  with  other  parameters,  the  influence  of  the  threshold 
energies  for  ionization  on  the  ionization  rates,  and, 
therefore,  should  probably  be  modified  when  the  threshold 
energies  for  ionization  are  explicitly  taken  into  account. 
Since  there  exists  no  theoretical  groundwork  to  indicate 
how  this  should  be  accomplished,  the  author  chose  not  to 
modify  the  functional  dependence  of  the  ionization  rates, 
but  merely  introduce  scaling  constants  c   and  c   so  that 


the  adjusted  ionization  rates  a   and  6   are  given  by 
a  =c,a  and  B  =c2^'  where  a  and  3  are  the  "reference", 
conventional  ionization  rates  given  by  van  Overstraeten 
and  de  Man  (1970) .   The  constants  c   and  c   are  adjusted 
to  give  the  best  possible  agreement  between  the  calculated 
and  the  experimentally  obtained  values  of  DC  multiplication. 

Since  the  adjusted  ionization  coefficients  would  be 
useless  if  c   and  c     were  determined  individually  for  each 
diode,  the  author  chose  to  adjust  c   and  c   on  diode  1N752, 
whose  DC  calculations  are  most  sensitive  to  the  values  of 
c,  and  c   because  they  need  to  attain  the  largest  values 
of  DC  multiplication  without  diverging.   Once  this  was 
done,  calculations  were  performed  on  the  remaining  diodes, 
using  their  respective  values  of  the  threshold  energies 
for  ionization,  as  determined  in  Chapter  IV.   The  results 
of  these  calculations  are  shown  in  Figures  47  to  49  under 
the  label  of  "calculated".   The  discrepancy  between  the 
calculated  and  the  experimental  values  of  DC  multiplication 
occurring  for  low  values  of  M  is  attributed  to  the  tri- 
angle approximation  for  the  junction  electric  field, 
which  gets  progressively  worse  with  decreasing  reverse 
voltage  since  the  diode  SCR  narrows  and  gets  progressively 
closer  to  the  changing  doping  densities  of  the  junction 
metallurgical  transition  region.   In  that  region  the 
triangle  approximation  predicts  stronger  electric  fields 
than  are  actually  present.   This,  in  turn,  leads  to  larger 


calculated  values  of  ionization  probabilities  and,  hence, 
larger  than  observed  values  of  DC  multiplication.   In 
spite  of  this  discrepancy,  the  DC  calculations  performed 
by  explicitly  taking  the  threshold  energies  for  ioniza- 
tion into  account  and  employing  the  adjusted  ionization 
rates  are  in  considerably  better  agreement  with  the 
experimental  data  than  are  the  results  obtained  by  con- 
ventional calculations. 

The  values  obtained  for  c   and  c   are  quite  reasonable 
also:   c -,=3,  c  =5.   In  fact,  it  was  on  the  basis  of  the 
magnitude  of  c„  that  the  author  dismissed  as  unlikely  the 
2.85±0.2  eV  value  for  the  hole  threshold  energy  for  ioniza- 
tion, for  when  this  value  was  employed  in  the  DC  calcula- 
tions it  was  not  possible  to  obtain  reasonable  agreement 
between  the  calculated  and  the  experimental  values  of  DC 
multiplication  for  values  of  c   as  large  as  twenty. 

It  was  also  observed  that  when  the  best  agreement 
between  the  theoretical  and  experimental  results  was 
obtained,  the  electron  and  hole  ionization  contributions 
to  the  DC  multiplication  decreased  monotonically  on 
successive  transits,  as  shown  in  Figure  43.   Since  the 
ratio  of  successive  terms  remained  essentially  constant 
regardless  of  which  terms  were  picked  (with  the  exception 
of  the  beginning  few)  this  implies  that  the  ionization 
probabilities  for  electrons  and  holes  were  essentially 


the  same.   This  information  is  used  in  Appendix  B  to 
simplify  the  derivation  of  the  expressions  for  the  magnitude 
of  the  low  frequency  spectral  density  of  noise  current 
produced  in  a  multiplication  process  resulting  from  no 
more  than  two  ionizations  per  carrier  transit  across  the 
diode  SCR. 


CHAPTER  VI 
CONCLUSIONS 

The  principal  result  of  this  study  is  a  detailed  under- 
standing of  the  mechanism  of  reverse  conduction  transition 
from  tunneling  to  impact  ionization  multiplication  in  low 
breakdown  voltage  silicon  p-n  junctions.   It  is  concluded 
that  the  transition  occurs  when  the  energy  of  the  tunneling 
electrons  exceeds  the  electron  threshold  energy  for  ionization. 

If  the  DC  multiplication  factor  is  redefined  by 
M  =  (I  /I  )-l,  and  the  functional  form  of  the  expression  for 
the  magnitude  of  the  noise  current  spectral  density  is  changed 
to  Si  =  2ql  [l+f(M)],  then  the  electron  and  hole  threshold 
energies  may  be  obtained  from  slope  transitions  in  the  plots 
of  log  f (M)  vs.  log  M.   The  electron  threshold  energy  for 
ionization  determined  in  this  manner  is  found  to  be  in  excel- 
lent agreement  with  the  value  of  1.8±0.1  eV  found  by  Moll 
and  van  Overstraeten  (1963) .   However,  the  hole  threshold 
energy  of  1.9±0.15  eV  disagrees  with  their  value  of  2.4±0.1  eV, 
Moreover,  it  is  found  that  the  threshold  energy  for  ionization 
is  not  constant,  but  increases  slightly  with  increasing  junc- 
tion width. 

The  noise  data  may  also  be  used  to  estimate  the  average 
minimum  distance  between  ionizing  collisions.   This  is 
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accomplished  by  obtaining  the  junction  width  at  the  reverse 

voltage  corresponding  to  the  onset  of  the  transition  from 

no  more  than  two  to  no  more  than  three  ionizing  collisions 

per  electron  or  hole  transit  across  the  diode  SCR,  and  dividing 

that  distance  by  three.   The  results  are: 

Diode  type   Average  min.  distance   Average  min.  distance 
for  electron  ionizing   for  hole  ionizing 
collisions  collisions 

1N749  192  A  199  A 

1N750  213  A  224  A 

1N751  223  A  234  A 

1N752  242  A  250  A 

When  the  threshold  energies  for  ionization,  the  total 
potential  at  the  onset  of  carrier  multiplication,  and  the 
value  of  (N  -N  )N  /N   are  interpreted  in  terms  of  the  transi- 
tion model  (Figure  35),  the  values  of  N  and  N   can  be 
determined.   The  doping  densities  on  both  sides  of  a  step 
junction  may  thus  be  evaluated  from  diode  noise  and  C-V  measure- 
ments.  This  technique  could  be  used  in  some  cases  for  the 
evaluation  of  impurity  solubilities  in  silicon. 

Conversely,  when  the  junction  doping  profile  and  the 
threshold  energies  for  ionization  are  known,  the  transition 
model  may  be  used  to  predict  the  reverse  potential  for  the 
onset  of  impact  ionization  multiplication.   This  furnishes 
a  method  for  distinguishing  the  tunneling  section  from  the 
multiplication  section  of  the  diode  reverse  V-I  characteristics. 


Besides  determining  the  onset  of  carrier  multiplication, 
the  threshold  energies  for  ionization  dominate  the  mechanism 
of  the  multiplication  process.   Because  they  allow  no  more 
than  two  or  three  ionizations  per  carrier  transit  across 
the  diode  SCR,  theyrender  Tager's  (1965)  and  Mclntyre's  (1966) 
theories  of  noise  in  avalanche  diodes  inapplicable  to  the 
diodes  of  this  study.   This  is  amply  demonstrated  by  the 
severe  discrepancies  between  their  theoretical  results  and 
the  experimental  data,  and  by  the  considerable  improvement 
in  agreement  between  the  data  and  the  theoretical  results 
obtained  when  an  upper  limit  of  two  ionizations  per  carrier 
transit  across  the  SCR  is  explicitly  imposed.   Indeed,  for 
low  multiplication  values,  a  noise  theory  based  on  the  explicit 
constraint  of  no  more  than  one  ionization  per  carrier  transit 
across  the  diode  SCR  gives  considerably  better  agreement  with 
the  experimental  data  than  do  the  theories  of  Tager  (1965) 
or  Mclntyre  (1966) . 

Moreover,  unless  the  effect  of  threshold  energies  for 
ionization  is  explicitly  taken  into  account,  a  finite  value 
of  DC  multiplication  is  calculated  when  no  multiplication 
is  experimentally  observed.   This  discrepancy  does  not  arise 
when  a  simulation  program  incorporating  the  threshold  energies 
for  ionization  is  used  for  calculating  the  DC  multiplication. 

The  many  details  of  the  transition  mechanism  thus  lead 
to  the  conclusion  that  the  conventional  treatment  of  carrier 
multiplication  phenomena,  based  on  the  equation 
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dJ  (x)    dJ  (x) 


a(x)  J  (x)+B(x)J  (x)       (6.1) 


dx       dx  n         p 

is  not  applicable  to  carrier  multiplication  in  low  break- 
down voltage  diodes. 


CHAPTER  VII 
RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  first  undertaking  in  further  work  on  this  topic 
should  be  the  removal  of  the  triangular  electric  field 
approximation  from  the  DC  multiplication  simulation  program. 
The  electric  field  should  be  obtained  from  the  Poisson 
equation  and  the  impurity  doping  profile.   If  this  correction 
does  not  remove  the  discrepancy  between  the  experimental  and 
the  calculated  values  of  DC  multiplication  at  low  reverse 
voltages,  the  cause  of  the  discrepancy  should  be  investigated. 

When  good  agreement  between  the  experimental  and  the 
calculated  values  of  DC  multiplication  is  obtained,  the  cal- 
culations should  be  extended  to  a  case  of  greater  practical 
importance:   the  reverse  breakdown  characteristics  of  base- 
emitter  junctions  of  bipolar  transistors.   Because  of  the 
very  heavy  doping  densities  encountered  in  base-emitter 
junctions,  it  will  probably  be  necessary  to  take  bandgap 
narrowing  into  account.   Moreover,  if  the  reverse  DC  V-I 
characteristics  of  the  base-emitter  junctions  or  of  the  low 
breakdown  voltage  diodes  are  to  be  calculated  directly  from 
impurity  doping  profiles,  it  will  also  be  necessary  to  refine 
the  theory  of  tunneling  breakdown  (Chynoweth   et  al. ,  1960; 
Tyagi,  1968). 
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Furthermore,  the  program  should  be  modified  to  calculate 
carrier  multiplication  resulting  from  a  single  pass  across 
the  space  charge  region.   This  would  permit  calculation  of 
carrier  multiplication  in  short  channel  FET's  and  other 
small  devices  in  which  one  carrier  type  is  somehow  removed 
so  that  it  does  not  participate  in  the  multiplication 
process. 

Since  the  threshold  energies  for  ionization  play  an 
important  role  in  the  simulation  program,  accurate  calcula- 
tions of  DC  multiplication  in  structures  of  arbitrary  lengths 
and  arbitrary  impurity  doping  profiles  will  require  accurate 
knowledge  of  the  threshold  energies  under  these  conditions. 
Consequently,  it  will  be  necessary  to  carry  out  additional 
DC  and  noise  studies  to  establish  means  of  predicting  the 
threshold  energies  for  ionization  in  such  structures. 

The  effect  of  different  junction  widths  and  different 
electric  field  profiles  could  be  investigated  further  with 
the  aid  of  the  recently  developed  replacement  series  for  the 
1N746  to  1N752  diodes  employed  in  this  study.   The  new  series 
of  diodes,  manufactured  by  alloying  aluminum  onto  the  lightly 
doped  p  layer  of  a  p-n   diode,  break  down  more  abruptly  than 
their  older  counterparts.   The  more  abrupt  breakdown  is 
obtained  from  the  resulting  p  -p-n   doping  profile  (author's 
conjecture!),  which,  for  any  given  reverse  voltage,  widens 
the  space  charge  region  and  lowers  the  peak  electric  field. 
The  lower  peak  electric  field  decreases  the  tunneling  current 


and  the  wider  junction  sets  up  a  more  favorable  environment 
for  impact  ionization.  Thus,  when  impact  ionization  multi- 
plication finally  sets  in,  it  occurs  more  abruptly  than  in 
the  conventional  p  -n  alloy  diodes  employed  in  this  study. 
Since  the  new  diodes  should,  for  a  given  value  of  DC  multi- 
plication, exhibit  more  noise  than  their  older  counterparts, 
this  conjecture  should  not  be  difficult  to  verify. 

Moreover,  it  would  be  possible  to  fabricate  narrower 
junctions  by  employing  stronger  p-type  material,  such  as 
boron  doped  aluminum  (Chynoweth   et  al.,  1960).   The  electric 
field  profiles  in  such  p   -n  diodes  would  be  more  asymmetrical, 
which  should  permit  impact  ionization  multiplication  at  lower 
reverse  voltages.   This  should  reduce  the  likelihood  of 
multiple  ionizations  per  carrier  transit  across  the  space 
charge  region,  and  make  the  noise  data  more  amenable  to 
interpretation  in  terms  of  the  already  developed  noise  theory. 

Due  to  phonon  involvement  in  the  carrier  transport  and 
energy  balance  mechanisms,  the  threshold  energies  for  ioniza- 
tion should  also  be  determined  as  function  of  temperature. 
Low  temperature  measurements  might  also  shed  some  light  on 
the  mechanism  of  phonon  collision  involvement  in  the  ioniza- 
tion and  multiplication  processes.   Tunnel  diode  DC  V-I 
characteristics  measured  at  4.2°K  revealed  structure  which 
reflected  the  cooperation  of  different  energy  phonons  in 
the  tunneling  process  (Chynoweth   et  al . ,  1962;  Logan  and 
Chynoweth,  1963;  Logan   et  al.,  1964).   Such  structure  might 
also  be  observed  in  the  DC  or  noise  data  of  these  diodes. 


In  conjunction  with  the  experimental  work,  the  noise 
theory  should  be  extended  to  the  case  of  a  multiple,  but 
finite  number  of  ionizing  collisions  per  carrier  transit 
across  the  diode  space  charge  region.   This  might  be  possible 
by  combining  the  author ' s  method  with  the  approach  employed 
by  van  Vliet  in  the  derivation  of  avalanche  photodetector 
noise  (van  Vliet,  1967) .   Incorporating  the  quantum  efficiency 
and  the  Boson  factor  in  the  primary  current  would  also  make 
the  theory  useful  for  evaluating  avalanche  photodetector  per- 
formance. 

Finally,  the  improved  multiple  ionization  noise  theory 
could  be  used  to  extract  from  the  noise  data  information  on 
the  ionization  probabilities.   This  would  be  of  considerable 
help  in  undertaking  a  theoretical  study  of  the  adjusted 
ionization  coefficients.   Such  a  study  might  also  lead  to 
new  approaches  to  the  treatment  of  other  aspects  of  carrier 
transport  in  small  devices  for  which  the  conventional, 
differential  equations  approach,  with  its  implicit  assumption 
of  unlimited  scattering,  is  not  valid. 


APPENDIX  A 

MULTIPLICATION  NOISE  RESULTING  FROM  NO  MORE  THAN  ONE 
IONIZATION  PER  CARRIER  TRANSIT  ACROSS  DIODE  SCR 


The  derivation  of  expressions  for  the  magnitude  of  the 
low  frequency  spectral  density  of  noise  current  produced  in 
impact  ionization  multiplication  resulting  from  no  more  than 
one  ionization  per  carrier  transit  across  the  diode  space 
charge  region  proceeds  as  follows.   As  indicated  in  Chapter  IV 
the  total  number  of  carriers,  N,  collected  at  the  diode  ter- 
minals in  one  second  is  given  by 


N  =   En. 
i-0  x 

where  n~  is  the  number  of  electrons  tunneling  in  one  second, 

and  n.,  i=l,  2,  3,  ...  is  the  number  of  hole-electron  pairs 

produced  in  one  second  on  the  i-th  carrier  transit  across  the 

diode  SCR. 

If  it  is  assumed  that  n   =an  ,  n  =bn~.  =abn  ,  n  =an  =a  bn~, 
etc. ,  where  the   denotes  averages  and  a  and  b  represent  the 

ionization  probabilities  of  electrons  and  holes,  respectively, 

then 

N  =  n0(l+a+ab+a2b+a2b2+. . .) 

and  the  DC  multiplication  factor,  M  =  N/nQ-l  becomes 

M  =  *<l+b) 
1-ab 
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The  low  frequency  spectrum  of  the  current  fluctuations 
can  be  obtained  from 

S   =  2var  N 

N 

where  S   denotes  the  magnitude  of  the  spectrum  of  N,  the 
number  of  events  occurring  per  unit  time  interval  (van  der 
Ziel,  1970).   Since  I   =  qN, 

2 

Sj   =  2q  var  N 

r 

2  ~~2  -2 
=  2q  (N  -N  ) 


2q2{(Z   n.)2-(  E  n  )  } 
i=0  X    i=0  1 


=  2q^{  Z   Z  (n.n.-n.n.) 


i=0j=0 


i  D   ID 


To  evaluate  n.n.,  it  is  necessary  to  have  explicit  expressions 

for  n.  and  n..   These  may  be  obtained  as  follows.   If  the  emission 
i      D 

of  tunneling  electrons  is  pictured  as  a  sequence  of  l's,  then 
the  terms  n.,  j=0,  1,  2,  ...  can  be  pictured  as 
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where  l.R  j=l,  2,  3,  ...,  k=l,  2,  3,  ...  nQ,  denotes  the 
occurrence  of  an  ionizing  collision  on  the  j-th  transit, 
initiated  by  the  k-th  tunneling  electron,  and  where  the  O's 
denote  the  absence  of  ionizing  collisions.   The  transitions 
from  1  to  0  occur  independently  and  at  random.   From  the  above 
tabulation,   it  is  apparent  that  n.  is  given  by 

nQ  ajr;  j  odd 

n.  =   £  a.  b.  a.  b,   ...  or 
D    r=1  lr  2r  3r  4r     ^ .  .  even 

where  a.,  and  b.,  ,  j=l,  2,  3,  ...  ,  k=l ,  2,  3,  ...  nQ  are 

1  or  0,  for  eventhough  a.,  and  b.,  are  allowed  to  fluctuate 

JK  .    3K 

independently  and  at  random,  II  a   b,   retains  the  two  essential 

k=l  Kr  Kr 
features  of  the  above  table,  namely,  that  (i)  if  all  preceding 

entries  in  the  k-th  column  are  1,  then  the  next  entry  can  be 

1  or  0,  and  (ii)  that  if  any  of  the  preceding  entries  is  0, 

then  the  next  succeeding  entry  must  be  0.   The  average 


n0  air   n0  ajp 

"l-J  ■  <r^alrb2r  ■•■«  "lalpb2p...    , 

ir   *  DP 

may  now  be  evaluated  by  ensemble  averaging.   Averaging  over  the 
m^  subensembles  containing  n   =  N  elements 
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Since  all  a's  and  b's  are  equal,  then  if  both  i  and  j  are  even 

n~TTN  =  N(aV    (b2)2    (a)       2       (b)    2 
j>i 

+N(N-l)(a)      2      (b) 
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Moreover,    since  a   and   b  are   1   or    0,    a   =a   and   b  =b.      Therefore, 

N     _4)_4)         .(*¥-) -i*& 

n±n.   =  Na  *    b    +  N(N-l)a      b      . 
j>i 
Finally,  averaging  over  all  subensembles 
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However,  since  the  tunneling  electrons  have  a  Poisson  distri- 
bution, for  which  var  n_  =  n_ 


-  -(2}-(2} 
(n.n. -n.n.)  =  nna    b     ;   i,j  even 
l  3   l  j  0 

j>i 


The  same  expression  is  obtained  if  j  is  even  and  i  is  odd. 
If  j  is  odd  and  i  is  odd  or  even, 

(n.n.-n.n..)  =  nQa     b  2 

If  i>.j ,  the  above  expressions  hold  provided  j  is  replaced  by  i 
Consequently , 
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Therefore, 
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If,  however,  the  probability  of  ionizing  collision  on  the 
first  transit,  a,,  is  not  equal  to  a,  the  probability  of 
ionizing  collision  on  succeeding  transits,  then 

(nTnT-n-.n.)  =  n^a^"1*  ;  j>i 

(n-n.-n.n.)  =  n^1^     ;  j<i 

and 

"*2  -2     

(nQ-n0)  =  nQ 


Consequently, 
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APPENDIX  B 

MULTIPLICATION  NOISE  RESULTING  FROM  NO  MORE  THAN  TWO 
IONIZATIONS  PER  CARRIER  TRANSIT  ACROSS  DIODE  SCR 


The  derivation  of  expressions  for  the  magnitude  of 
the  low  frequency  spectral  density  of  noise  current  pro- 
duced in  impact  ionization  multiplication  resulting  from 
no  more  than  two  ionizations  per  carrier  transit  across 
the  diode  space  charge  region  proceeds  as  follows. 

It  is  assumed  that  carrier  multiplication  taking  place 
during  one  transit  across  the  diode  SCR  is,  in  general, 
of  the  form 


*   II- 

1 

\ 
\ 

al 

\ 
1              \ 

»1*2 

[ 

*2 

*2 

/ 

/ 

1 

where  "1"  denotes  the  "primary"  carrier  initiating  the 
multiplication  process,  "x"  denotes  an  ionizing  collision, 
a,  and  a_  denote  the  probabilities  of  the  first  and  second 
ionizing  collision,  and  the  dashed  lines  indicate  carriers 
of  opposite  charge  to  the  "primary"  carriers. 

If  the  above  diagram  is  extended  to  include  several 
transits  across  the  diode  SCR,  it  takes  on  the  appearance 


of  a  "tree"  from  which  it  becomes  apparent  that  the  multi- 
plication process  can  be  decomposed  into  an  infinite  number 
of  "branches"  of  the  form 

-2-2      -3-2 


na,      na.a.    nana_    na.a.     na.a. 
x Ux  L-£x    1   \   x   1   2.  x L*.  ...  etc, 


which  is  precisely  the  up  to  one  ionization  per  carrier 
transit  across  the  SCR  process  for  the  case  of  unequal 
probabilities  for  holes  and  electrons.   (It  should  be 
observed,  however,  that  in  this  analysis,  the  ionization 
probabilities  for  electrons  and  holes  are  assumed  to  be 
the  same.   The  evidence  supporting  this  assumption  is  dis- 
cussed in  Chapter  V.   Even  if  this  assumption  were  not 
strictly  valid,  however,  Figure  40  indicates  that  the  results 
obtained  for  unequal  ionization  probabilities  do  not  differ 
substantially  from  the  result  obtained  for  equal  ionization 
probabilities  in  the  one  ionization  per  carrier  transit  case. 
And  since  the  up  to  two  ionizations  per  carrier  transit 
across  the  SCR  case  is  a  composition  of  the  one  ionization 
per  carrier  transit  cases,  the  final  result  should  not  be 
very  sensitive  to  the  difference  in  ionization  probabilities 
for  holes  and  electrons . ) 

To  reduce  the  number  of  "branches"  that  need  to  be 

considered,  it  will  also  be  assumed  that  a_<<a  ,  so  that 

-   -i-k 

branches  beginning  with  a.,a|a2,  k  >_2 ,  will  be  neglected. 

This  assumption  is  justified  in  Appendix  C  .   Moreover, 


7  5 


since  the  probability  of  ionizing  collision  is  smaller  on 
the  first  transit  than  it  is  on  succeeding  transits,  the 
probability  of  the  first  ionizing  collision  on  the  first 
transit  across  the  diode  SCR  will  be  denoted  by  a,  . 
According  to  the  conduction  model  of  Figure  35,  the  prob- 
ability of  the  second  ionizing  collision  on  the  first  transit, 
a   ,  for  the  range  of  voltages  where  the  up  to  two  ionizations 
per  transit  across  the  SCR  theory  may  be  considered  valid, 
is  zero.   This  follows  from  equation  3.13,    with  1.8±0.1 
replaced  by  2V  .   The  results  of  this  calculation  are 

Diode  type  Reverse  voltage  for  which  a?->>0 

1N749  4.89-5.06 

1N750  4.85-5.00 

1N751  4.90-5.06 

1N752  4.95-5.10 

The  reverse  voltages  given  above  are  greater  than  the  voltages 
corresponding  to  the  first  slope  transition  in  Figures  36  to 
39,  thus  confirming  the  contention  that  for  the  range  of 
voltages  where  the  up  to  two  ionizing  collisions  per  carrier 
transit  theory  is  valid,  a   =0. 

Under  the  condition  stipulated,  the  "tree"  for  this 
multiplication  process  is  given  in  Figure  50.   The  heavy 
arrows  denote  electron  ionizations,  the  light  arrows  denote 
hole  ionizations,  and  the  I kj|  denote  "branch"  numbers.   To 
avoid  excessive  clutter,  only  the  beginnings  of  each  "branch" 
are  indicated.   Thus,  "branch"  I 111  is  the  branch  given  by 


aHaia2 


allaia2 


-  -4-2 
ailaia2 


etc. 

All  branches,  therefore,  continue  vertically  downward, 
and  their  continuation  is  indicated  by  dashed  arrows.   The 
beginning  of  each  new  branch  is  indicated  by  a  horizontal 
translation  to  the  right,  followed  by  the  downward  arrows 
of  that  branch. 

The  analysis  of  this  multiplication  process  proceeds 
as  follows.   Let  N  be  the  total  number  of  carriers  collected 
at  the  diode  terminals  as  a  result  of  multiplication  initiated 
by  nQ  tunneling  electrons.   Then, 

N  =  Im. 
i  x 

where  m.  denotes  the  number  of  carriers  produced  by  the  i-th 

branch.   Hence, 


var  N  =  I  var  m .+211  (m .m . -m .  m.) 
i         lj  J 

j>i 


But 


m.  =   E  n. , 
1   k=0  lk 

where  n.,  is  the  number  of  carriers  produced  in  the  i-th 

branch,  on  the  k-th  transit  across  the  SCR.   Following  the 

procedure  of  Appendix  A,  n.,  is,  in  turn,  represented  by 

IK 


n.,  =   E    (product  of  a,,,  a  '  s  and  a  's) 

where  a....,  a  '  s  and  a  's  are  1  or  0,  and  the  product  function 
of  a.,,  a  ' s,  and  a  's  reflects  the  past  history  of  the  r-th 
primary  carrier,  in  the  i-th  branch,  on  the  k-th  transit 
across  the  diode  SCR.   Thus  the  problem  of  up  to  two  ionizing 
collisions  per  carrier  transit  across  the  SCR  is  reduced 
to  a  sum  of  an  infinite  number  of  up  to  one  ionization  per 
carrier  transit  analyses. 

To  obtain  the  expression  for  M,  the  branch  contributions 
are  summed  as  follows: 


N  =  Em.  -m1(J  +  l^    +  _EQmojo 


+   E    Em...+   Em._+   E    Em.. 
i=2  j-1  0l^    j=2  3°    i=2  j-1  ^ 


resulting  in 


n  ri  +  i   f  1+^  )  ■  *H5?52(1+52) 

no[1  +  an(  -  -  } 


1-ala2     (l-a1a  ) (1-a  ) 


5ll*l52(1+52)    ,  5llala2(1+52) 


(l-i1I2)  (l-i1)2    (l-i1a2)  (l-a;L)2 


ill5l(1+i2)       illili2(1+i2)   , 
d-a1a2)  (1-a^    (l-51a2)  (1-a^)2 


N 
Hence,  M  =  —  -  1  becomes 


a,a~    1+a,      a,  (l-a,+2a,  a_)  (l+a0) 

m  =  a    [(1+-M.)  (  _  _  )  +  — — _1 2 1  2  _ — ] 

l-ax   1_a1a2       d-a1)  (l-a1a2) 
To  calculate  the  results  for  var  N,  the  procedure  of 
Appendix  A  is  duplicated  to  obtain 

-  -   -2- 
a.. (3+5a1+a1a--a.a«) 

var  m  =  1  +  -±± ± — ±-S — ±-£ — 

1  (l-a^r 

and  the  variances  of  the  branches  b,  and  b2  whose  averages 
are  given  by 

bx:   l+a1+a1a2+a12a2+a12a22+. . . 
and 

b2:   i1+£1i2+i12i2+a12i22+... 


They  are,  respectively 


_  _   _  2- 
l+3a,  +3a1a:)+a,  a~ 

var  bx  = — _  _  2   — 


and 


a,  (l+3a2+3a1a2+a;.a2  ) 
var  b9  =  = 

Observing  that  the  variances  of  branches  whose  averages 
are  represented  by 

illi1^2k(1+5l+5lVili2+  ■••> 


-   -  j~  k.-   -  -   -  2-   -2-2 
allal   2  *al   la2   1  a2   1  a2  +  *•• 


are  given,  respectively,  by 


a,  ,  a -,  a2   var  b. 


ini1Da2k  var  b2 


and  repeating  the  summation  performed  to  obtain  N  =  Em. 

i  x 


Z  var  m .  =  1  + 
i 


all(3+5al+ala2-ala2) 
(1"ala2)2 


+  allala2r1+3al+3ala2+ala2 
+  _ [ __ 


(1_ala2) 


allal d-S1+25lS2)  r1+3a2+3ala2+ala2\ 
|_ j 


(1-a^ 


(l-a1a2) 


It,  therefore,  remains  to  compute  2Z  E  (m.m.-m.m. ). 

i  j     J 
j>i 

Fortunately,  it  is  necessary  to  compute  only  seven  (m .  m  -m~. m . )  '  s , 
since  the  remaining  can  be  obtained  by  inspection.   The  pro- 
cedure employed  in  calculating  these  is  illustrated  for  the 


case  of  (m10m11-m10  m^) 


<m10mll"m10mll)  =  -^  .f1(n10inllj"n10inllj) 


where 


and 


j  terms 
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Here,  ^_  denotes  identically  the  same  factor  in  n, Q.  and 
nlli  '  anc^  the  f actors  bracketed  by        are  to  be  treated 
as  one,  since  they  correspond  to  i=l  or  j=l.   Since  this 
approach  is  identical  to  that  employed  in  Appendix  A,  it 
may  be  observed  there  that 


^lO^llj^lO^llj5  =  (n10inllj)|r=s 
since  n02-n02=n0  . 

~~2~    ~2  "~ 2  

Moreover,  when  r=s,  a  . =a, , ,  a,=a,,  and  a  =a2.   Consequently, 


(m10mll"m10mll)  =  ±f1  .^^loAlj^: 


a11a1a2(l+a2)  (5+a2-4a1a2) 

n, 


d-i152)2 


The  other  six  (m.m . -m .m . ) ' s  of  importance  are 


(m11m12-m11m12)  =  a^a^i—J^-)2 

1_aia2 


a11a1a2(l+a1)    a2(l+a]L) 

{m10m020-mlOm020)    =         n    -   -    , [3+n    -   -    ,] 

(l-a1a2)  (l-a.^) 


,        l+a0  (1+a.) 

2      _} 


(m10m02rmllm021)    =   allala2(1    -    -    }  [2+77- 

1-a.,  a_  (1- 


-    alV 
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^ll^lS2(1+al)  (1+a2} 


(mlim020-miim020)  " 

(l-a-j^) 


3r  ,  1+a2  ,  „  .  1+a2 


(m020m021-m020in021)  =  allala2  (-^-^)  (1+7^^-) 

1"aia2     1'ala2 


-4-  ,  1+a2  ,,„.  1+a2 


(m2dn2i"m2on,2i)  =  aiiaia2(^-^-)(2+r-^-) 

1"ala2     1"aia2 


All  other  (m.m  .-m.m. )  •  s  may  be  obtained  from  these  by  observing 


that,  for  example, 


(mlimij-miimij)  =  al(mliml(j-l)-mliml(j-l)) 
i>l,  j>3 


Generating  additional  (m.m. -m.m. ) • s  in  this  manner  and  adding 
the  results 
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E        E    (m,  .mn  . -m.  ,mn  .  ) 
i=0  j-1     Xl  ^     Xl  1] 
j>i 


n0anaia2  (1+a2)  (5+a2-4a1a2) 
(l-i1)  (l-i1a2)Z 


nQallala2,    1+a2    ,2 
2_( — — — -) 

(l-a]L)z   l-axa2 


.^0(m10m02j-m10m02j)  = 


n0allala2(1+a2)(3+a2~2ala2) 
(l-a1)  (l-i1a2)2 


n0allala2 ^1+al^  (3_2aia2+a2^ 

(l-i152)2 


-  -   -3-2    -      - 
n  an  na,a0  (1+a, )  (l+a_) 


.Zn  .Zn(mlim02j_mlim02j) 
1=1  j=0  J 


l0°HT2,i,ul'  %  _2' 


(l-a1)  (l-a1a2) 


+   H0allala2(    1+a2    }2 
(1-i-L)  2       l"aia2 


I    (m,  ^m-  .-m,  nm0  . )    - 


HOaliai (1+a2)  (3+a2"2aia2) 


i=0 


10m2j    Ul10lll2j 


(1-a^  d-a1a2) 


(l-a1+a1a2) 

Z zr~z — 2 

(l-a^  (l-a1a2) 


°°        °° x         n0allala2(1~al+ala2)  ,      1+a2N2 

E         E     (m      m      -m      m      )    =   1 5 ( — it-' 

11    2^       lx    2=>  (l-a1)2  1"ala2 


i=l    j=0 


Wj-0      °21    °2:      °21    °2:  tt-S.) <!-«.=-> 2 

j>i  L  12 


.    n0aHaia2f    1+a2       ,2 
(l-a1)         l-^1a2 


-    -      -4_ 

™       °° nnanaia9      1+a9  1+aT 

Z        Z    (m2im2i-m2iS2i)    =  -  < — )(2+ 2-) 

i=0    j=1      2!    2D       2i    2D  (  -    -  x_- 

j  >i  i  ±    2  1    2 


n0aliaia2/    1+a2    ,2 

= — T~ ( — =""=-) 

(l-a1)         1_a1a2 


jf0(m20m02j-m20m02j) 


"0allala2(1+al) (1+a2) 


(l-a1a2)2 


+   K0aHaia2(      1+i2)2 
(l-a1)         l-axa2 


*        *     (rn2im02i-ra?imn^)    -   "0aHala2(1+al)(1+a2) 

i=l    i  =  o  D        x    u^  = ^~3 — 5 

J  (l-ax)  (l-a^r 


"0ailai6aL    1+a2    ,2 

= — 2~( — -   -    i 

(l-a1)         1-aia2 
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00 -    -      _2-2      1+ai 

.!n(m020m03j~m020m03j)    =   n0allala2(1     -    -    ) 
J-U  12 


_    _      -4-2         -  - 

n0allala2 (1+al)  (1+a2) 

(l-i1)  (l-i1i2)2 


E         £     (m_„  .m_ -j  • -rom  ,ni-n-3  •  )    = 


Vll5li2(1+il)(1+52) 


.^   ^ -02±«03 j   «02i-03j'    -       (1.5i)(1.5i52)2 


Vllala2,    1+*2    ,2 

z — 2~~ ( — -   -    > 

(l-a1)         1_aia2 


.Zn(m020m3j_m020m3j) 
3  =  0  J  J 


"0allala2(1+al)  (1+a2)  (1_a1+aia2) 
(1-a^)  (l-a^)  2 


00    °° n0allala2(1_al+ala2) 

i!l  jf0(m02im3j-m02ira3j)  -  ^J-l 


1+a         2 
•  (      _         ) 

1"ala2 


N: 


jf0(m20m03j_m2  0m03j 


n0aliaia2(1+ai)(1+a2) 


(l-a1a2) 


+  H0allala2(  1+a2  j 2 
V2 


"0illili2(1+il)(1+i2: 


Z         E     ^2im03j-in2im03j)    =  (1_-    -    ,2 

i=l    j=0  1    2 


,505115152(    1+i2    }2 

d-S1)       i-i1i2 


H0ill5l(1-5l+il52)/       1+^2,2 


P:  .E  <m2om3i-m2om3j)  = ttt"; (7TT) 

3  =  0  J  (l-a^  l-a1a2 

*Vll5l*2(1~*l+*l*2) 
Q:  E  E     (m2im3n_m2im3i)    =  -       2 

i=i  j=0    2l  3D    2l  3:  d-a  r 


<^->2 


Finally, 


2E  E (m.m.-m.m. )  =  2 [A+ — —  (B+C+D+E+F+G+H+I) 
ii   1]   1]         1-ii 
j>i 

+  - T(J+K+L+M+N+0+P+Q)  ] 

(1-^) 

After  a77,  aT  ,  and  a^  are  calculated  using  the  DC  multiplica- 
tion program  of  Appendix  C,  the  above  results  are  evaluated, 
and  (var  N)  -  1  is  plotted  vs.  M  in  Figures  44  to  46. 


APPENDIX  C 
DC  MULTIPLICATION  SIMULATION  PROGRAM 

The  DC  multiplication  simulation  program  performs 
calculations  of  carrier  multiplication  resulting  from 
impact  ionization  in  p-n  step  junctions.   Its  input  variables 
are: 
NSETS    :   number  of  diodes  whose  DC  V-I  characteristics 

will  be  calculated 
NVREV    :   number  of  reverse  voltages  for  which  the  DC 

multiplication  will  be  calculated 


XNA 

XND 

ALPHA1 

Bl 

ALPHA2 

B2 

ALPHA3 
B3 

PHI 
DV 


N  -N   =  net  doping  density  on  p-side  (cm   ) 
N_  =  net  doping  density  on  n-side  (cm   ) 
electron  ionization  rate  constants:   adjusted 
ionization  rate  is  given  by  a  =ALPHAl  e  E 
hole  ionization  rate  constants:   adjusted  ioni- 
zation rate  for  E  <  4  x  10   V/cm  is  given  by 

_B2   ~ 
B,  =  ALPHA2  e  F~ 

hole  ionization  rate  constants:   adjusted  ioni- 
zation rate  for  E  >  4  x  10   V/cm  is  given  by 

B3 
B2  =  ALPHA3  e  E 

<J),  junction  built-in  potential 

Xj+Ax. 
dV  =  -q  /        E(x)  dx,  where  Ax   is  the  length 
k 

of  the  k-th  interval  in  the  E (x)  vs.  x  diagram 
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TV 
EV 
HV 
SLOPE 


tunneling  potential 

electron  threshold  potential  for  ionization 

hole  threshold  potential  for  ionization 

the  tunneling  current,  I  ,  is  given  by 
SLOPE 

I, 


t  -  i//V 


VREV     :   reverse  voltage  for  which  DC  multiplication 

will  be  calculated 

The  different  portions  for  the  program  perform  the 

following  computations.   Statements  29  to  66  divide  the  p 

and  n  sides  of  the  electric  field  triangle  into  equipotential 

intervals,  DXN(K),  DXP (K) ,  and  compute  the  electric  field 

values,  EN(K),  EP (K) ,  at  the  end  points  of  the  intervals. 

Statements  67  to  95  number  these  quantities  consecutively, 

starting  on  the  p  side  with  the  first  interval  and  the 

first  electric  field  value.   Statements  96  to  164  compute 

the  approximate  values  of  the  ionization  probabilities 

for  electrons  and  holes,  dp  (k)  =  XINTN(K)  and 

n 

dP  (k)  =  XINTP(K),  in  each  interval,  DX(K).   The  spatial 
distribution  of  the  tunneling  current,  dl  (k)  =  TI (K) ,  is 
computed  in  statements  173  to  219. 

The  electron  initiated  multiplication  is  calculated 

in  statements  220  to  241,  251  to  266,  and  276  to  287. 

TV 
Statements  220  to  225  shift  the  tunneling  current  by  ^ 

intervals  to  obtain  the  free  electron  current.   (In  "con- 
ventional calculations"  TV  is  set  to  one  DV. )   The  book- 
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keeping  algorithms  which  trace  the  electrons  and  holes 
back  and  forth  across  the  SCR  and  compute  the  DC  multi- 
plication on  each  transit  are  contained  in  statements  237 
to  241  and  262  to  266.   (If  the  available  storage  is 
sufficiently  large,  the  encremental  increase  in  current, 
DIE(N)  and  DIH(N),  resulting  from  carrier  multiplication 
in  each  interval  can  be  replaced  by  DI(K,N),  where  K  is 
the  transit  number  and  N  is  the  interval  number.   Combining 

£  DI(K,N)  with  the  interval  number  as  function  of  distance 

k 

then  yields  the  distribution  of  carrier  multiplication  as 

function  of  distance  across  the  SCR.) 

Statements  242  to  250  and  267  to  275  retrieve  and  save 
the  electron  multiplication  in  each  interval  on  the  10-th 
and  12-th  transits,  and  the  hole  multiplication  on  the  11-th 
and  13-th  transits.   These  quantities  are  used  later  in  the 
program  to  compute  the  values  of  a^a-,  and  a~,  the  first, 
second,  and  third  ionization  probabilities.   (There  is 
nothing  special  about  these  transits.   Any  other  corresponding 
set  of  transits  for  which  the  multiplication  has  achieved 
a  steady  state  condition  would  do  as  well.) 

The  multiplication  initiated  by  holes  left  in  the 
valence  band  when  electrons  tunnel  to  the  conduction  band 
is  calculated  in  statements  293  to  338.   They  are  essentially 
mirror  images  of  the  statements  handling  the  electron  initiated 
multiplication.   The  total  DC  multiplication  is  then  the 
sum  of  the  electron  initiated  and  the  hole  initiated  multi- 


plication.   It  is  denoted  by  SUM  and  computed  in  statement 
342.   This  ends  the  DC  multiplication  simulation. 

The  remainder  of  the  program  is  devoted  to  the  calcula- 
tion of  the  ionization  probabilities,  a,,  a_  ,  iL ,  and  to 
the  evaluation  of  the  final  results  of  the  up  to  two 
ionization  per  SCR  transit  noise  theory. 

The  electron  ionization  probabilities  are  computed 
in  statements  345  to  421.   The  probability  of  the  first 
ionization,  a,,  is  computed  by  taking  the  multiplication 
due  to  holes  on  the  11-th  and  13-th  transits  and  calculating 
the  electron  ionizations  on  the  12-th  and  14-th  transits, 
allowing  no  more  than  one  ionization.   This  is  performed 
in  statements  345  to  356  and  379  to  390.   The  square  root 
of  the  ratio  of  the  first  electron  ionizations  on  the  14-th 
transit,  El,  to  the  first  electron  ionizations  on  the  12-th 
transit,  Dl,  is  a,.   If  this  calculation  is  repeated  starting 
with  the  first  electron  ionizations  on  the  12-th  and  14-th 
transits  (the  results  of  the  computation  just  completed) , 
as  done  in  statements  358  to  367  and  391  to  401,  a„  is 
obtained.   Repeating  the  calculations  once  more  in  state- 
ments 368  to  378  and  402  to  412  leads  to  a^. 

The  corresponding  hole  ionization  probabilities  are 
computed  in  the  same  manner,  but  starting  with  the  electron 
multiplication  on  10-th  and  12-th  transits.   These  calcula- 
tions are  carried  out  in  statements  452  to  519. 
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Because  the  two  ionizations  noise  theory  does  not 
distinguish  between  electron  and  hole  ionization  proba- 
bilities, it  is  now  necessary  to  evaluate  the  "average" 
ionization  probabilities  to  be  used  in  the  noise  theory. 
They  are  simply  the  square  root  of  the  product  of  the 
electron  and  hole  ionization  probabilities,  as  given  in 
statements  558  to  560. 

According  to  Appendix  B,  the  probability  of  the  second 
ionization  on  the  first  transit  is  zero  for  all  reverse 
voltages  for  which  the  probability  of  the  third  ionization 
is  zero.   Hence,  a,-,  =  S(2),  the  probability  of  ionization 
on  the  first  transit. 

Now,  that  the  ionization  probabilities  a..  ,  ,  a,  and  a~ 
are  known,  the  theoretical  expressions  for  noise  resulting 
from  up  to  two  ionizations  per  carrier  transit  across  the 
SCR  may  be  evaluated.   This  is  performed  in  statements 
561  to  577.   This  completes  the  program. 

The  calculations  described  above  yield  the  following 
ionization  probabilities  for  the  1N751  diode: 
Diode  reverse  voltage         a,  ,         a-,         a~ 
5.05  0.485       0.909      0.114 

5.0  0.474       0.897      0.107 

4.95  0.462       0.886      0.0994 

4.9  0.451       0.874      0.0926 

4.8  0.428       0.851      0.0796 

4.7  0.403       0.829      0.0676 


Diode  reverse 

voltage 

all 

al 

a2 

4.6 

0.378 

0.807 

0.0567 

4.5 

0.353 

0.786 

0.0469 

4.4 

0.327 

0.766 

0.0381 

4.3 

0.299 

0.746 

0.0304 

4.2 

0.273 

0.728 

0.0237 

4.1 

0.245 

0.710 

0.0180 

4.0 

0.218 

0.693 

0.0133 

3.9 

0.190 

0.677 

0.0094 

3.8 

0.163 

0.661 

0.0064 

The  range  and  functional  dependence  of  the  above  values 
are  typical  of  the  other  diodes  employed  in  this  study. 
Consequently,  the  approximation  that  a„  <<  a, ,  invoked 
in  the  derivation  of  the  up  to  two  ionizations  per  SCR 
transit  noise  theory,  is  justified. 
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DIMENSION      XN(600),VNf600),EMAX(600),DXN(600),EN(600),DXP(600), 
1XPC60O), VP(600).EP(fcOO),fc(    900),XInTn(    900}»XlNTP(    900), DX(    900), 
2VREV(305.TIC900). 
3DIE(1200),OTHC1?00) 

REAR    13,    CiSPTS,NVB£\ 

FORMATC2I?) 


3(30!),SOIN(900),SOI(900),SSOIC900), 


DO  1000  NS=1,  ^SETS 

RFAOlO,  xua.xm5.ai PHA1,R1, ALPHAS, R2,ALPHA3,B3 

10  FnRMAT(2F9.?,6E7.3) 

READ  11,  PHl,DV,TV,EV,HV#SLOPE 

11  FORMATES. 3, F10.7) 

REAP    12.     fvPEV(l5,I=l,NvREV) 

12  FORMAT(16F5,3) 


PRINT    lU»(vSEV(I),Iel,NVPFV2 
FORMATC'IVREV1 


IH  FORMATC1VREV'       /,14(1X,F5.3,2X)) 

EPSPFLsll.7 

EPsr=B.Pfe*-u 

G=l  .6E-1 9 

wl=?rjBT((2.*EP3PEL*EPSO/C)*C(XNA*XN0)/(XNA*XND))) 

DP  1000  NVsi,mVrev 

Wrwl*SQPTfPMi+VREV(UV)) 

Ga2.*S|_PPE/Vl 

E«AXX=C2.*SnRT(PHI  +  v/REV{tW)))/Wl 

IF(PMAxx.r,E.  1.2EA)  GO  TO  98 

PDT.NT  97 

97  FORMATC1X, 'MAXIMUM  FIELO  VALUE  TOO  SMALL1) 
STOP 

98  XU\=tXUA/(XUA  +  X'oD))*W 

C      COMPUTE  OX  AND  E  FOR  N  SIDE 
XN(1)=XNN 
^(1)  =  V»  ' 
■MAXCIK. 
(3VNK/PV 
IF  (K-1J  101,99,99 

99  DO_100  fl=l,K 

DXnHi)  =  Xm(II)*(i  -SQOTCl.-DV/VUCII))) 
EK'(Tl5=PMAXriI)*(f  ,-DXN(lf)/XM(II)) 
Vurjj)=VN(lT)-DV 
XU?JJ)=VN(IT)-0XU(II) 
Em1x(JJ)=FN(H) 

100  CONTINUE 
DXNtK+n=2.*(VNW-K*DV)/ENCK) 
GO  TO  102 

101  CONTINUE 
JJsl 
DXNf JJ)=2.*VNN/EMAXX 

C  COMPUTE    OX    AND    F    FOR    P    SIDE 

102  XPC1)=XPP 
IS  VP(1)=VPP 
19  EMAX(1 )s£MAXX 
iO  M»VPP/DV 
il  IF     (M-l)     201,199,199 
i2                      199  DO    200    1=1, m 
S3  Jel>l 

Dxpfn=xp(ii*(i.-sQRT(i,-oy/vP(i))) 

EP(I)sEMAxtl)*(!.-DxPnI/XP(I)) 

VP(J)=VPCIJ-0V 

xpm=xpM-oxpm 

EMix(J)=EP(l) 
200    CONTINUE 
GO  TO  202 

DXP(J)=2,*VPP/EMAXX 
GO  TO  203 

202  CONTINUE 
DXPCm+1)=?.*CVPP-m*Dv)/EPCM) 

203  CONTINUE 
C      RENUMBER  OX 


201 


IaM+1 

O  300  m=1,T] 
x(n)=dxpJ>j3 
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0071  JN  =  .!N-1 

007?  300    CHNTIMUF 
0073  JJJ»M*K. 

007U  Ji 

0075  3 


0076 
0077 


\W2 


DO  310  NsJA.JJJ 
OXtNj.DXNfjt) 

JI=JI*1 

0079  310    CONTINUE 

„„.,.  C      RFNUMBER  E 

888?  319    5^=5H   J  ^1,32,,319 

0082  josh 

0083  5n  320  NEsl.jQ 
008g              E(nE)=EPCJN$ 
0085              JNsJn-1 

8S|^       us  mwu** 

oin  322   ]£.&   J  "°'"0.322 

0090  JEEbM+K+1 

0091  jp-i 

2221  on  330  WE=JE,JEE 

0093  E(NE)»EN(jPr 

009<i  .rp  =  .TP*i 

0095  33 


OF  APPROXIMATING 


?      £2M5"II  IHNTZATTDN  SURlNTEGRALS  FOR  ELECTRONS 

I  §ubintegr,al*t  AK'D  3T0P  p6ints  for  tRaPez°i°  Method 

B«AM*x1(B1,p2,B3) 

0097  IIcM*l 

0098  I!J:km 

0099  33j    DO  3un  1=1,11 

0100  KESTAR=I 

0  01  BPE=B/E(I) 

8iHS  3*0    ^N??NUFE-100-}  G°    T°  "l 

OJ  oa  341    DO  3a?  1  =  1,  TIJ 

0105  KFSTOPsJJ.t-l 

O|0*  BOE=R/E(KE3T0P) 

0|07  IF(RQE.LE.l00.5  GO  TO  3u3 

0108  3«2    CONTINUE 

0j09         3«3    ItIkESTArI?  *ND  XINTP'3  ™  ZER0  FnR  FIELD  VALUES  WHICH  ARE  TOO  SMALL 
OllJ  IF(Il".LF.O)  GO  TO  3U5 

0111  KK=M*1 

DO  3«U  1  =  1, KK 

345         II=KESTOP+2 

IFUJJ-KESTOP-l    .LE.O)    GO    TO    347 

OP~3-*>    1  =  1, KK 

xintn(it)=5 
xintp(ii)=o 

3«6  k2J}2cix-lt-01  go  to  3u7 

S,.-,         TRIANGLE    APPROXTMATION    FOR    FIRST    AND    LAST    NON    2FRO    XINTN. 
347  XlNTN(KESTAR)  =  (r.X(KESTAR5*ALPHAl*FXPC-Bl/E(KESTAR)n/? 

?^INf^STr-1P+"  =  (nXCKF3TOP+n*ALPMAi;EXP(-H1/F    KESTOP)J)/2. 
C  I2AP^0TO    METHfin    II3E0    FOR     REMAINING       XINJN'S  ' 

IF(KESTnp-KESTAP    .uT.    1)    GO    TO    U01 

Ii=*esTpp-i 

DO  «00  I=KEST4R,II 
J=I*1 
«00    ^^^^  =  0xfJ)*4LPH*l*(EXPt-Bl/E(n)*EXP(-Bl/E(J)))/2. 

UP 


2  Cn«PiiTF«    FIRST    novi    ZERO    XINTF 

01         IFCEfKFSTARjjLE,    u.FSl    GO    TO    a02 

HA3« 
XINTP  (KESTAR)  =  (OX rKESTAR)*ALPHA2*EXPt-B2/E(KESTAR)))/2. 


xJ,N'TPfKFSTAP5  =  CnxrKESTAR)*ALPHA3*EXP(-B3/E(KE8TAR)))/2< 
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C 
403 


404 

I 

412 


450 
451 


504 
505 


506 
507 


mi 


COMPUTE  LAST  NO 

IF(F(*ESTOP+l)     .IF.    4.E5)    GO    TO    404 

XINTP(KESTOP+ns(PXCKFST&P+l)*ALPHA3*EXP(-a3/E(KESTfiP5n/?. 

GO  TO  412 

XTNTPCKF3T0P*l)=(0XCKFSTUP*n*ALPHA2*EXP(-B2/'E(KEST0P))J/2. 

USE  TRAPEZOID  APPROXIMATION  FDR  REMAINING  XINTP'S  WHILE  CONSIDERING  FIELD 

VALUES. 

IFfKFSTnp  .  KESTAP  .LT.  1)  GO  TO  501 

Dn"50b  IsKESTAR.H 

IF(Efl)  .LE,  4.E55  GO  TO  451 
XINTP(J)=DX(J)*ALPH43*CEXP(-B3/ECin*EXP(-a3/F(J)))/2, 

:(J)«ALPHA2*(EXP(-B2/E(n)+EXP(-B2/ECJ)))/2. 

CONTINUE 

CONTINUE 

PRINT  19,  VREVfVI 

FOH"At('1',t2,'vrfv  =  ',F5.3) 

PRINT  15,  XNA,XVD,ALPMA1,B1,ALPHA2,B2(ALPMA3,B3 

format (/T4, «xna',ti<5, 'xm6' ,T32, ' alpha ^ ,juo,1het a «,/4(ix, eio, u f4X)/ 

lT32.2(E10.4.5X2/T32.2fEl0.4,5X)) 

PRINT  16,  PHI,0V,TV,EV,HV 

FORMATC/1X, 'PHI  =  ',F5.3,5X,'DV  x  •  , FS .3 , 5X, • TV  *    » ,F5.3,5X, • Ev  » 
1».F5.3,5X,'HV  s  I.F5.3) 

PRIMf  18,  M 

FORMAT(/,lX. «M  =  ',131 

PRI>'T  17,  F"AXX 

FORMAT?/, IX, 'PEAK  FIELD  IN  JUNCTION  a    •  ,E10 ,4,3X, ' VOLTS/CMI ) 

FSrmItT/Iix, 'JUNCTION  WIDTH  s  • ,E10,4,3X, 'CM'3 

PRINT  STATEMENT  For  DX, XlNTN, XINTP, AND  E  GOES  HERE 

NTO=TV/DV+l . 

NTT=mTO*1 

NFO=EV/DV+l. 

NET=NEO*l 

NHOBHV/DV+1. 

JHn  =  J.!J  +  NHO 

JMO=JJJ-l 

JMHPsJJJ-NHO 

COMPUTE  DISTRIBUTION  OF 

COMPUTE  START  AnD  STOP  P 

H3«li 

DO  502  T=1,TI 

KESTAR=I  . 

BOEsG/ECI) 

IF  (BOE.LE.50.)  GO  TO  503 

CONTINUE 

DQ  504  1=1, IIJ 

TOP) 
IF  (BOE.LE.50.)  Go  to  505 

CONTINUE 

SET  TUNNELING  PROBABILITY  TO  ZERO  FOR  SMALL  FIELD  VALUES 

TTsKESTAR-1 

lFCII.LE.0D  GO  TO  507 

DO~506  1  =  1, KK 

iU  nuee-°5  go  to  5o7 

CONTINUE 

II=*ESTOP*2 

IF  (JMO-KESTOP.LE.0)  GO  TO  5C 

KKrK+1 

00  50$  1=1, *K 

TI(in  =  0. 

11=11*1 

IF  (JJJ-II.LT.O)  GO  TO  509 

CONTINiJF 

Continue 

TRlANO.iF  APPROXIMATION  FfR  FIRST  AND  LAST  TUNNELING  SUBInTEGRAl 

TTfKFSTAR)  =  {DX(^f.STARi*FXP(-G/ECKESTARjn/2. 

TI(i»ESTnp*l)  =  (DX(KESTOP*n*EXP(-G/E(KEST6P)J)/2. 


TUNNELING  PROBABILITY 

0INTS  FOR  TRAPEZOID  METHOD  OF  APPRQX  SUBINTEGRAL3 


JUF 
TUNNELING  PROBABILITY  PRINT  STATEMEN1 
ELECTRON  INITIATED  MULTIPLICATION 
SHIFT  TUNNELING  PROBABILITY  TO  OBTAIf 


I  JIJ^!kA^G.P,322*!lLlTY  PRINT_STATEMENTS  GO  HERF 


o|o2       c     if*^^jf8o*KpEp^KR:p?.?rgoN?sG5?riNTEGB*Ls 

0|0T  ^ijo^KESTAR,!! 

SliO        510    J^J^X(J»*(EyP(-G/eCI))*EXP(-6/E(J))))/2. 

0211  511     CONTINUF 

02l?  C      siJHTpboZE  T'-INNEL1*'G  P0Q"ABILITY  FOR  EACH  3UP1NTERVAL 

o?n  OH  51?  1=1  ,JJJ 

QiU  „    SllMTP  =  SI'MTP*TI(I) 

OflS  512    CONTINUE 

0216  IIsKESTOP*l 

2?}2  DO  513  I=KESTAR,IT 

22J2  ^  ,    TICT)  =  TT(n/8UMT- 

0219  513    CPNTINi-- 
UNNEL3 

_lectrc 

0220  C      g^^^^L^C'PBOBAPiLifrTO'OBTAlN  TUNNELING  CURRENT 

022?  550    CONTINUE 

22"  DO  551  NSMTT,JJ.T 

0**g  __.    DTEfK'5  sTICN-NTO) 

Of t^  551    CONTINUE 

2|2§  Sdisn. 

Hll  Dn    552    N3NTT.JJJ 

0230         C      STOM^M  "ITI*TE0  MULTIPLICATION 
02.11  kkbk 

0*2?  KMOSK-1 

0233  OP  560  N=1,NE0 

0?-3u  DlHCN)sO. 

0235  560    CONTINUE 

0235  Hi=o, 

0|37  Dn  5Sl  N=MET,JJJ 

2?4S  Xl=BI+0IEfN-NEC)+2.*DIH(N-NE0) 

0239  OIHCN)  =XTNTNCN)*XI 

2?*?  RlaXI-OIh(N) 

0?ai  561    CHMTINUF 

H~,$  „,     lp  (K-10)  563,520,522 

2|«?  520    DO  52t  Nsl,Jj5    ' 

gijfS         521    CONTINUE  l  ' 
02<i6  GO  TP  563 

IF  fK-l?)  563,523,563 

On  5?u  nsi,jjj 

VP(N)bDTHCN) 
524    CONTINUE 
563    SfKlsO. 

OP  56?  NsNET.JJJ 

SCK)=.SCK)+DTHCN) 
562    cnNTTMUF 

IPC^CO.LE.l.E-B)  GO  TO  581 

KPn=K+l 

KKaKPO 

R?^72  n°Jmmo,jjj 

DIE(N)c0. 
570   CnNTINUE 
Rlafl. 

On  571  N=NHn,jMn 
XlaPI+DTH(JHG-N)+?,*DTE(JHO-N) 
=  Xf.NTPCJJJ-N)*XI 


III 


DIEMJJ-N) 
RIaXI-DTEtJJJ-M) 
571    CnNTINUE 

IF^-101  586,587,589 

587  On  58*  K'=1  ,JJJ 
XNCN)=DTE(N) 

588  CPNTfMuF 
GP  TO  SP6 

I§?    JPCK-12)  586,590,586 
590    On  591  ^»1,JJJ 

XP(N)=DIt(N) 


97 


0275  591         CONTINUE 

037*.  586         S(KPP)=0. 

0277  DP    S72    msi.jmho 

0278  S(KPC1sS(KPD)+0TEfN5 

0279  572    CONTINUE 

0280  IF(S(KPO)/Sf2).LT.0.002)GO  TO  581 

0281  580    CONTINUE 
0??2  581    CONTINUE 
02e3  8TEST=Sf2) 
0?*u  Al 1  =  S  f 25 


SNsO. 

DO  585  NS2, 
SwsSN+SCN) 

CONTINUE 

FnRMlT(/>/i:iX. 'ELECTRON  INITIATED  MULTIPLICATION' 
"10.1) 


U?«6 
0287 
028*         585 


0289  PRINT  21,  3^ 

0290  21     FORMAT(///,l 
t'SN  =  ' .ElO.u) 

0291  POU'T  20,  (*■,  Sfn,  N=i.KK) 


Sfn,  Nsi.K 

0292         20     FnRMAT(//(T?/'Sf '.13,')  =  '.E10.4)) 
C      HOLE  INITIATED  "ULTIPI ICA"-- 


.  .  TATtO  "ULTIPt  ICATION 
TRANSFER  TUNNELING  PROBABILITY  TO  DlH(N) 

0293  DO  600  N31,.!JJ 

0294  DjHfN)  iTII«) 

0295  600    CONTINUE 

0296  sm  =  o. 

0297  DO  601  nsi,,jjj 

0298  S(nsS(l)*OTH(N) 

0299  601    CONTINUE 
COMPUTE  HOLE  INITIATED  MULTIPLICATION 
DO  MO  K»2, 100,2 


8^8£ 

0302 


KMP=K-1 

0303  DO  610  NsjMHOfJJJ 

0304  DIE(n)b0. 

0305  610    CONTINUE 
0^06  RI=0, 

0307  DO  611  nxwho,Jmo 

0  308  XtsRI+niH(JMO-N1+2.*DTE(JHO-N) 

0309  DTE(JJJ-N)  cXINTP(JJJ-N)*X1 

0310  RI=XI-DIEf JJJ-N) 

0311  611    COMTIVUF 

0312  S(K1=0. 

0313  DO  612  Nal,.JMHO 

0314  SCtOsSCO+DTECN) 

0315  612    CONTINUE 

0316  IF(9(K).LE.l.E-5)  GO  TO  641 

0317  KP0=K*1 

0318  KKsKPO 

0319  DO  620  Net ,wEP 

0320  0IH(N)=0. 

0321  620    CnNTINUF 

0322  RI=0, 

0323  DO  621  NsNET,JJJ 

032"  XT  =  RT  +  DIE(M-NEOU2.*DTH(K-NEO) 

0325  DIH(N)    SXINTNCN^XI 

0326  RI=XI-0IHCN) 

0327  621    CONTINUE 

0328  3(KPOl  =  fl. 

0329  DO  622  nsmft.JJ.T 

0330  sckpq)=?(kpo)+dth(n) 

0331  622   Continue 

0332  IFfSjKPOJ/STEST.LT.O.Ol)  GO  TO  641 

0333  640    CONTINUE 

0334  641    CONTINUE 

0335  8P=0. 

0336  DO  650  M»2,KK 

0337  SP=SP+S(N) 

0338  650    CONTINUE 

0340  24           FPRmIt!///??X,'holE    INITIATED    MULTIPLICATION' ,//,T2, ' 8P    *    »,E10.4) 

03'U  PRINT    20,     (N,S(N),                          N=1,KK) 

03<i?  SIImsSN  +  SP 

cas  print  ?*,   sum 

03a4  26           F0fiMAT(/////,T2,'SUM   a    SN    ♦    SP    »    ',El0.4) 


JJ>5  DO    720    N=1,JJJ 

0346  DIE(M)sXMfN5 
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0347  DTHtNian. 

0348  720    CONTINUE 
0JU9  721    HT=0. 

035?  DO730  wawET.JJJ 

035?  XIsPI+OTEfN-NEO) 

03sa  ri=xi-oth(nj 

0355  Dleni+OTHfNj 

0357  Dn  731  n=1,JJJ 

0358  DIE(N)=0. 

0359  731    CONTINUE 

0360  RT=P. 

°0\U  0n=720  m=mET,JJJ 

o?63  8l5?It2IWS2^91, 


0364 


0419 
0«|0 
0421 
0422 


DIE(lO  =  XlNTN(N)*XI 


0365  RIcVI-OIE(N) 

§366  D?=C2*OIEfM) 

0367  740    CONTINUE 

036«  DO  74}  Nsl,jJJ 

0369  DIH(n5=0. 

0370  741    CONTINUE 

0371  glfO. 

83?|  0n~750  N=NET,JJJ 

0376  RIsXI-DIH(N) 

0377  D3sD3*OTH(N) 

0379  00  760  Nsl.JJJ 

§380  OIE(n)=XpH) 

0381  DIH(N)=0. 

0382  760    CONTINUE 

0383  RJsO. 

j^gg  EoB770  N=NET,JJJ 

a\ll  XT=Ri+niEfN-NEO) 

§?§?  6iM(N)SXlNTM(N)*XI 

0388  RlrXI-OTMfN) 

0389  ElsPl+OIHtMj 

0390  770    CONTINUE 

0391  00  771  nsi,JJJ 
039?  DIEfMjsO. 

0393  771    CONTINUE 

0394  RIsO. 

0395  E2S2*. 

0396  00  780 
0397 
0398 

0400  E2=E2+DTEC 

0401  780    COUTINUF  m 
quo?  DO  781  Nsl.JJJ 
04f>3  DlH(O  =  0. 

0404  781    CONTINUE 

0405  ?I=2* 

gSSf  50=7*0    M.NET.JJJ 

0408  XI  =  "I  +  DIE(N-NECI5 

5a59  DIH(Ni=XINTN(N)*) 

0410  RIcXI-DIHtNJ 

0412  790    CONTINUE 

0413  RATlni=n2/Dl 

0414  RATI0?=03/D1 
0(415  R*TI03  =  E2/E1 


0n  780  n=net,JJJ 
XIsRI*DIH(n-nE05 

DIE(N)=XINTN(N)*XI 


R»TI04=E3/E1  , 
A1=SQRTCE1/D1) 


A?soiTIOJ 
R2sA? 
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04?3  789    FnRMAT(//,i*, 'noise  PANAMETE3S  COMPUTED  F30*  ELECTRON  IONIZATION') 

0424  PRINT  799, All 

0425  799    FOWMATC/,  J  X,'AHs  ',E10.4)      

0«2*  PRINT  801,01 

0427  800    FORMATC/,lX,'0l  s  ',£10.45 

0U28  PRINT  801,02 

0«29  801    F0«*AT(/,5 <. '02  =  ',E10.*n 

0430  Point  «n?.r>3 

0431  802    FnwATC/,lx.'D3  =  '»E10.4) 
0'«3?  PRI^T  803, Fl 

C-133  803    FDR"ATC/,1X, 'El  =  ',E10.4) 

0434  PRINT  804. t? 

0435  804    FORMAT(/,lX,'E2  a  ',E10.4) 

0436  PRINT  805. E3 

0437  805    FORMAT(/,f X, 'E3  =  ',£10,4) 


044  0 

0441  807 


0438  PRINT  806,RATlni 

0439  806    FHRMATC/,1 X, '02/01  =  ',E10.4) 
POIHT  8"7, RATIO? 
FPR»AT(/,f X. 'Dj/01  =  ',E10.4) 

044?  PRINT  808,*ATIO3 

044?  808    FORMAT?/, IX, 'E2/E1  =  ',E10.4) 

04(1<I  PRINT  809, RATI  04 

04«5  809    FnRMATC/,lX, 'E3/E1  =  ',E10.4) 

0446  PRINT  810, A1 

04U7  810    FORMAT(/,lX, 'Al  s  ',E10.4) 

04UB  PRINT  811, A? 

04U9  811     FPRMAT(/, 1 X, ' A2  s  ',E10.4) 

0450  PPI»T  812, RATIOS 

0451  812     FHRMAT(/,1  X,  '  AU/A1=  ',E10,4) 

83H  B9E?3Wi!WJJ 

0454  0IHfN)=0, 

0455  830    CONTINUE 

0456  RI=0. 

0457  OlsO. 

0458  On  831  M=NHn,jMn 

0459  XT=RI+0IEtJH0-N) 

0  4*0  DThrjJJ-^)=xiNTPC.7JJ-k,)*XI 

0461  RI  =  XI-0I,1CJ.TJ-N} 

046?  QlsCI+ClHf JJJ-N) 

0463  831    CONTINUE 

0464  DO  832  N=1,.TJJ 
0U65  DTE{N)=0. 

0466  832    Continue 

0467  RT=0. 

0468  O?=0. 

0469  DO  840  NSMHO.JMn 
0470 
0471 
047? 
0473 
0474  840     CONTINUF 

0"i75  DO  841  N=1,.IJJ 

0476  DTHfNlcO. 

0477  841    CONTINUE 

0478  RI=0. 
0"79  03  =  0. 

0480  Do  850  n=nho,Jho 

0481  XI=PI+PTECJwn-N) 
DTH(JJJ-N)  =  VINTP(.IJJ-W)*XI 
RI  =  XI-0IHCJ  TJ-N1 


0482 
0483 


XIxPT*DTMf JM3-N) 
DIF(JJJ-N)  =  XlNTPt,TJJ-N)#XI 
RI  =  YI-DIECJ.TJ-N) 
Q2s02  +  DIE(J.TJ-N) 


0484 
0485 

04S6  DO  flfcO   K=1,JJJ 

0487  OlfcfN)sVP(N) 

04fl«  0IHfM)s0. 

0489  860    CONTINUE 

0490  RI=0. 
0491 
0492 

0493  XlsRI*DIE{JHC-N)         ww 

0  4RU  DIH(JJJ-N1=XIMTP(JJJ-N)*XI 

0495  RT=XT-OIHf  J.tJ-MJ 

0'196  TlsTlfOTHf  J.TJ-N1 

0497  870    £!2^Tr|L-,F 

0498  00  871  n=i,JJJ 


0099 

0500 

871 

0501 

0502 

0505 

050a 

0505 

0506 

0507 

0509 

880 

0509 

0510 

0511 
0512 

881 

0513 

0510 

0515 

0516 

0517 

0518 
0519 

890 

0520 

0521 

891 

052? 

0553 

900 

052O 

0525 

901 

8ii7 

902 

0528 

0529 

903 

0530 

SHI 

900 

0533 

905 

0530 

0535 
0536 

906 

0537 

0538 

0539 

907 

05J0 

054] 

0542 

908 

0503 

0500 

0505 

909 

0506 

0507 

0508 

910 

0509 
0550 

0551 

911 

055? 

0553 

055O 

912 

0555 

0556 

913 

0557 

0C58 

0559 

0560 

0561 

056? 

0563 

0560 

0565 

0566 

0*67 

0568 

0569 

057(1 

ml 

0573 

imv»& 


RU 

T?  =  0. 

On  880  NSNHD.jMn 

XTsPl+DlHfJMO-N) 

DTEfJJJ-M=*INTP(,TJ,T-M)«XI 

RIsXT-OTEfJJJ-NJ 

T?sT2*DIt(JJJ-N) 

CONTINUE 

On  Pel  K'st  ,jjj 

DtHCN)sO, 

CONTINUE 

RT  =  0. 

13  =  0. 

OP  890  vsMMn.jMn 

xi=»i*oiErjHo-Ni 

OlHf  JJJ-M=YIMTP(.JJ,T-M)*XI 

RTsXI-DIM(J.TJ-N) 

T3  =  T3  +  DH(JJJ-Ni 
CONTINUE 

Fni^lTC//, IX, 'NOISE  PARAMETERS  COMPUTED  FROM  HOLE  IONIZATION') 

PRINT  900.01 

FOPMATC/,1 X, 'CI  =  »,E10.4) 

PRINT  901,0? 

Fn«"AT(/,lx, '02  =  »,Ei0.O) 

PRINT  9(12,03 

FOHmaTC/.IX, '03  s  '.E10.O) 

PRINT    9«3, T1 

FHRmatC/.IX.  'Tl    =    '.E10.O) 

PPPT    900,  T? 

FOR"AT(/,  IX,  'T2    =     ■  ,E10.O) 

PRINT    905, T3 

FPPMAT(/,lx,'T3    =    ',E10,«) 

RAT1=C2/01 

PRINT    906.     RATI 

FO«MAT(/,tX,  '"2/01     =     ',E10.0D 

RAT?=Q3/01 

PRINT  907.RAT? 

FnRMAT(/,!x, »  C3/01  =  ',ElO.O) 

RAT3=T?/T1 

PRINT  908.RAT3 

FPRVATC/,1X,  'T2/T1  =  ',E10.O) 

RATUST3/T1 

PRINT  909.RATO 

FTR^ATC/.IX, 'T3/T1     -     ',F10,U) 
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Figure  1.   Typical  voltage-current  characteristics  of 
commercially  available  diodes 
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Figure  22.   Reverse  V-I  characteristic:   diode  1N746 
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Figure  23.   Reverse  V-I  characteristic:   diode  1N747 
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Figure  24.   Reverse  V-I  characteristic:   diode  1N748 
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Figure  25.   Reverse  V-I  characteristic:   diode  1N749 


loDi \ 1 1 1 1 1 1 1 r 


iol 


10- 


5  102 


10 


O 


10' 


10" 


rl 


1N750 


J I I I 1 L 


0.4>4  0.43  0.52  0.56 

OL+1.01)"* 


0.60 


Figure  26.   Reverse  V-I  characteristic:   diode  1N750 
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Figure  27.   Reverse  V-I  characteristic:   diode  1N751 
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Figure  28.   Reverse  V-I  characteristic:   diode  1N752 
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Figure    29, 


Comparison  of  Tager's  theory  with  data 
diode  1N749 
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Figure  30.   Comparison  of  Tager's  theory  with  data: 
diode  1N750 
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Figure  31.   Comparison  of  Tager • s  theory  with  data: 
diode  1N751 
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Figure  32.   Comparison  of  Tager's  theory  with  data: 
diode  1N752 
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Figure    33 
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Figure  34.   Mclntyre's  theory 
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Figure    35.       Multiplication   threshold   model 
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Figure  36.   Multiplication  noise:   diode  1N749 
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Figure  37.   Multiplication  noise:   diode  1N750 
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Figure  39.   Multiplication  noise:   diode  1N752 
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Figure  40.   Multiplication  noise  resulting  from  no  more  than 
one  ionization  per  carrier  transit  "..j.oss  SCR 
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Figure  41.   Comparison  of  one  ionization  per  SCR  transit 
noise  theory  with  data:   diode  1N7^' 
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Figure  43.   Contribution  to  multiplication  on  n-th  carrier 
transit  across  the  diode  SCR 
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Figure  44.   Comparison  of  noise  theories  with  data: 
diode  1N749 
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Figure  45.   Comparison  of  noise  theories  with  data: 
diode  1N750 
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